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We discuss the QCD cvohition equations governing the high energy 
behavior of scattering amphtudes at the leading logarithmic level. This 
hierarchy of equations accommodates normal BFKL dynamics, Pomeron 
mergings and Pomeron splittings. Pomeron loops are built in the course of 
evolution and the scattering amplitudes satisfy the unitarity bound. 
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1. Introduction 

Over the last three decades, one of the main active fields of research 
within QCD has been the study of its behavior in the high energy hmit. In 
general, a scattering process will be considered as a high energy one, when 
the (square of the) total energy s of the colliding (hadronic) objects is much 
larger than the momentum transfer between them. At the same time 
one hopes to approach the problem via analytical methods, since in this 
limit there is the possibility of a large kinematical window s ^ ^ ^qcd 
where one can apply weak coupling methods. 

The first approach to the problem was done in the mid seventies when 
the BFKL (Balitsky, Fadin, Kuraev, Lipatov) equation [1], one of the cen- 
tral equations governing the approach to high energy, was derived. It was 
understood that certain Feynman diagrams in perturbation theory are en- 
hanced by logarithms of the energy and therefore they have to be resummed. 
This equation was solved in a special case (for the forward amplitude), and 
a total cross section growing as a power of the energy emerged. This growth 
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may not be so surprising (at least) a posteriori, since at high energies the 
wavefunction of a hadron can contain a large number of partons, mainly 
gluons, due to the available phase space for virtual fluctuations and due to 
the triple-gluon coupling in QCD. In general, QCD at high energy will be 
characterized by high densities and increasing (but constraint) cross sec- 
tions. 

Afterwards, in the early to mid eighties, it was realized that one needs 
to find a mechanism to tame the too steep increase of the partonic densi- 
ties as predicted by the BFKL equation and the concept of saturation was 
introduced as a dual description of unitarity in high energy scattering. The 
gluon density at a given momentum should never exceed a value of order 
0{l/as) and, as the mechanism to fulfill this saturation bound, a non-linear 
term was proposed [2] to be added to the (linear) BFKL equation. Later 
on a proof of that equation in a special limit (the double logarithmic limit, 
which nowadays is known not to be a good approximation for a high den- 
sity system) was given [3]. At the same time another step of progress was 
made, as the solution of the BFKL equation for non-forward scattering, or 
equivalently at fixed impact parameter, was obtained [4]. 

In the mid nineties, one can say that there was a major breaJsthrough. 
It was also the time when various subficlds were created, as different ap- 
proaches to address and approach the high energy problem started to de- 
velop. The color dipole picture [5-8] was formulated as a description for 
the wavefunction of an energetic hadron, a picture which goes beyond the 
BFKL equation in the sense that it also contains transitions between dif- 
ferent number of Pomerons in the multi-color limit (with a Pomeron de- 
fined, more or less, as the object which evolves with energy according to 
the BFKL equation). This allowed one to calculate higher moments of the 
gluon densities while at the same time the approach to unitarity limits could 
be studied. A program to calculate the vertices for these Pomeron transi- 
tions (beyond the large- A'^c limit) was started in [9]. Roughly at the same 
time a somewhat different problem was studied, namely the saturation of 
densities in the wavefunction of a fast moving large nucleus, due to the 
strong coherent classical fields generated by the large number of its valence 
partons [10]. Even though there was no QCD evolution in that model, the 
ideas introduced proved of great significance for what followed in the forth- 
coming years. Moreover, that period faced the first attempts to attack the 
high energy problem by the method of effective actions [11-13]. Finally a 
Hamiltonian, equivalent to an integrable system, was given for a particular 
configuration of n "reggeized" gluons [14,15]. 

Even more progress was seen during the last years of the previous decade 
and the beginning of the current one. The idea introduced earlier in [10] 
that the wavefunction of an energetic hadron can be described in terms of 
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strong classical color fields, was used properly in order to derive an equa- 
tion describing the evolution to higher and higher energies. This functional 
equation, called the JIMWLK (Jalilian-Marian, lancu, McLerran, Weigert, 
Leonidov, Kovner) equation [16-19], gives the evolution of the probability to 
find a given configuration of a color field associated with the hadronic wave- 
function. In general, the physical system which is supposed to be described 
by such an equation, like a fast moving proton, nucleus or a quarkonium, 
was called a Color Glass Condensate (CGC) [17], for reasons to be explained 
later. In turn, the JIMWLK equation can be used to derive equations for 
the scattering amplitudes of given projectiles off the evolved hadron, and 
the outcome is what is known by today as the Balitsky hierarchy [20] . This 
is a set of non- linear equations, which under a mean field and/or large- 
Nc approximation collapses to a single one [21] giving the evolution of the 
amplitude for a color dipole to scatter off the CGC. As a byproduct, the 
odderon problem (scattering with an odd number of gluon exchanges) in- 
troduced in the early eighties [22-24] was reformulated and extended to its 
non-linear version [25,26]. During that period another important accom- 
plishment was the calculation of the next to leading order correction to the 
BFKL equation, which was finally completed in [27, 28], while there was an 
effort to calculate the vertices for Pomeron transitions, their properties and 
consequences [29-33] , a task which is still ongoing. 

Until the beginning of the last year it was thought, at least in some part 
of the "high energy community" , that the JIMWLK equation was a more 
or less "complete" and self-consistent description of the high energy limit of 
QCD (at the leading logarithmic level). However, a calculation done within 
the dipole picture resulted in some large corrections [34] for the saturation 
momentum (the momentum scale at a given energy at which gluonic modes 
saturate). This result was clearly different from what was known up to 
that time, and the discrepancy was initially attributed to the difference 
between the JIMWLK equation and its mean field version. Shortly it was 
understood that these corrections are an effect of the low density fluctuations 
in the dilute, high-momentum, tail of the wavefunction [35] and in fact the 
significance of these fluctuations in the evolution of the system had been also 
observed and realized much earlier from numerical simulations within the 
dipole picture [36, 37]. The understanding that the JIMWLK equation does 
not describe properly the low density region of the hadronic wavefunction 
came soon [38, 39], and a generalization of the Balitsky equations at large- 
Nc was given. This new hierarchy does not allow any kind of mean field 
approximations and contains loops of Pomerons, in contrast to the Balitsky- 
JIMWLK one which contains only "one-way" Pomeron transitions. And in 
fact, the possibility to allow for the formation of Pomeron loops in the 
course of evolution is of crucial importance for a self-consistent approach to 
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unitarity. Not surprisingly, within this QCD description and under certain 
logical approximations, the results in [34, 35] were reproduced. Triggered by 
these observations, facts and derivations, there has been an ongoing effort to 
the direction of constructing a Hamiltonian formulation and a generalization 
of the equations at finite- A^c [40-45]. 

In these lectures we will start by introducing in the next section the 
BFKL Pomeron and the BFKL equation in coordinate space as was derived 
in the dipole picture. In Sec. 3 we will discuss its pathologies and in Sec. 4 
the concepts of saturation and Pomeron mergings will be proposed as the 
resolution to the problem, while at the same time a non-linear equation 
will naturally emerge. In Sec. 5 the Color Glass Condensate, the JIMWLK 
equation and the Balitsky equations will be presented, and in Sec. 6 the cal- 
culation for the energy dependence of the saturation momentum within this 
formulation will be given. In Sees. 7 and 8 we will try to explain which are 
the main problems encountered within the JIMWLK evolution and which 
is the most crucial missing element. In Sec. 9 the dipole picture will be re- 
viewed and the significance of fluctuations will be stressed, while in Sees. 10 
and 11 the new hierarchy will be derived, and then the way Pomerons loops 
are generated will be quite obvious. In the next three sections we will dis- 
cuss the approach to formulate the problem at the Hamiltonian level, a 
remarkable (possible) property of the Hamiltonian and the approach to the 
generalization of the theory beyond the multi-color limit. In the last section 
the saturation momentum will be revisited, since its energy dependence will 
be influenced by the low density behavior of the effective theory. As it is 
clear the presentation will be mostly based on the dipole picture and the 
Color Glass Condensate formulation of the high energy problem, and over 
the years there have been nice lectures and reviews for both approaches [46- 
56] . We shall not discuss at all any phenomenological aspects of the theory 
and the possible importance of high density phenomena in the relevant ex- 
periments. Nevertheless, we need to say that there have been successful 
qualitative and sometimes quantitative descriptions of the the low-x data 
in deep inelastic lepton hadron scattering (DIS) [57-60] and of the observed 
particle spectra in heavy- ion collisions [56,61-67]. 



2. The BFKL Pomeron and the BFKL Equation 

We start these lectures by introducing the concept and the significance of 
the BFKL Pomeron and by giving a heuristic, but also intuitive, derivation 
of the BFKL equation. Imagine that we want to measure the gluon distribu- 
tion of a generic hadron. One way to do this, is by probing the hadron with 
a small, in size, color dipole. If the dipole has a size r = \x — y\ ^qcd' 
where x and y are the coordinates of its quark and antiquark legs respec- 
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tively, it will probe the gluonic components of the hadron with momenta 
such that Q ~ 1/r ^ Aqcd- Clearly, the condition that the dipole be small, 
is dictated by the requirement that the QCD coupling constant be small and 
thus the problem can be approached by analytical methods. Such a probe 
is not only a good theoretical object due to its overall color neutrality, but 
it can also be "created" as a fluctuation of the virtual photon in deep in- 
elastic lepton-hadron scattering (DIS). Therefore, if one is able to calculate 
the dipole-hadron scattering amplitude, one can obtain the cross section for 
DIS, since the probability for the creation of the dipole is determined by 
a lowest order calculation in QED. The latter can be easily performed [68] 
and therefore we shall restrict ourselves to the analysis of the dipole-hadron 
scattering. 

Let us assume that the hadron, to which we shall frequently refer as 
the target, is right-moving, while the projectile dipole is left moving. We 
shall also assume that the target is energetic enough, so that a partonic 
description of its wavefunction is meaningful, while the dipole is slow enough 
so that it is "bare"; that is, there are no additional components, through 
higher order radiative corrections, in the dipolar wavefunction. 

At lowest order in QCD perturbation theory the dipole and the hadron 
will interact via a two-gluon exchange as shown in the left part of Fig. 1. 
Now let the hadron be boosted at very high energy and let be the (light- 
cone) longitudinal momentum of a "valence" parton of the hadronic wave- 
function^. Then, and as we show in Appendix A in more detail, the probabil- 
ity to emit a soft gluon with longitudinal momentum in the interval from k'^ 
to k'^ + dA;+, with /c+ = xp^ ^ 1, is proportional to agAk'^/k'^ = agdx/x. 
When X is small enough and due to the QCD triple-gluon coupling, the 
gluon with momentum fraction x can be generated through the intermedi- 
ate emission of one or more gluons which have their longitudinal momenta 
strongly ordered, i.e. they satisfy x <C a;„ ^ • • • ^ xi ^ 1, where Xi is 
the momentum fraction of the z-th gluon in this sequence of emissions of n 
intermediate gluons. This process, when compared to the direct emission of 
the soft gluon at A;+ from the valence parton at p"*", is of order^ 



where the logarithms have been clearly generated from the integration over 

^ Since we are interested in the high energy hmit, any possible masses will be always 
taken to be zero. 

^ We shall always let the coupling be fixed (only in Sec. 6 we shall briefly discuss some 
extensions to running coupling) and for the reasons explained earlier we shall consider 
it to be small. 
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Fig. 1. From lowest order in pQCD to the BFKL Pomeron. 



the available longitudinal phase space for the intermediate emissions. At 
high energies, the rapidity Y, defined as Y = ln(l/x), can compensate the 
smallness of the coupling Ug, so that OsY > 1. Thus one needs to resum 
all these (agY)^ enhanced terms and such a resummation gives rise to the 
so-called BFKL Pomeron. It is obvious in Eq. (1), that this procedure will 
lead to an exponential, in Y, increase of the gluon density and this will 
be verified shortly and in more detail when we also take into account the 
degrees of freedom in the transverse phase space. The interaction of an 
energetic hadronic target with a color dipole is shown on the right part of 
Fig. 1, where the summation over the gluon ladder represents the BFKL 
Pomeron or equivalently the small- a; components of the hadronic wavefunc- 
tion. Notice that this figure corresponds already to the square of a diagram 
in perturbation theory, since we are interested in determining the probabil- 
ity to find a given mode inside the hadronic wavefunction. 

Instead of performing the resummation of diagrams, one can equiva- 
lently study the evolution of the target wavefunction or of the scattering 
amplitude under a step Ay in rapidity. Furthermore, one can also view the 
last soft gluon as being emitted from the projectile dipole, which is a simple 
object and therefore its evolution can be easily studied. In the large-A'^c 
limit the soft gluon can be represented by a quark-antiquark pair, and thus 
the final system is composed of two dipoles {x,z) and {z,y), with z the 
transverse coordinate of the emitted soft gluon. As we show in Appendix A 
the differential probability for this splitting process is [5] 



where = agNc/ir, with the number of colors. In the above equation, 
dY = dln(l/x) = —dk^/k~^ represents the differential enhancement in the 
longitudinal phase space as shown in detail in Appendix A, while the kernel 
M xyz contains all the QCD dynamics in the transverse phase space. Now, 
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Fig. 2. The BFKL equation for the dipole-hadron scattering amphtude. 



each of the two final dipoles can scatter off the target and therefore the 
evolution equation for the imaginary part of the scattering amplitude T 
reads 



dY 



BFKL 27r 



j Mxyz [{Txz) + {Tzy) - {Txy)] = as^EFKL® (T^y) ,(3) 



where the last term is the virtual contribution arising from the normaliza- 
tion of the dipole wavefunction, and with the average to be taken over the 
target wavefunction, even though this is irrelevant for the moment. This 
is the BFKL equation [1] in coordinate space [5, 69] and its diagrammatic 
representation is shown in Fig. 2. It is free of any divergencies, since the 
potential singularities arising from the poles of the dipole kernel cancel. For 
example, when z = x the last two terms in the square bracket cancel each 
other, while the first term vanishes due to color transparency. In Appendix 
C we give a more rigorous derivation of Eq. (3), by using the corresponding 
BFKL equation for the density of dipoles which is derived in Sec. 9. 

The BFKL equation is a linear one, and therefore the solution can be 
obtained by studying the corresponding eigenvalue problem. Let us first 
look at the simplified case where the scattering amplitude is integrated over 
the impact parameter b = {x + y)/2 to give the total cross section, and 
is invariant under rotations in r. Then it depends only on the magnitude 
r = \x — y\ of the projectile dipole. As we show in Appendix D the solution 
to this eigenvalue problem is 

/CBFKL55r2(i-7)=x(7)r2^'-^\ (4) 

where 

x(7) = 2V(l)-V'(7)-V'(l-7), (5) 

with V'(7) the logarithmic derivative of the F-function. The eigenvalue func- 
tion x{l) is plotted in Fig. 3. Thus, one can write the solution to the BFKL 
equation as the superposition of the evolved with rapidity eigenfunctions, 
namely 

(T(r, y)) = ^ 1 ^ To(7) exp [asx{l)Y + (1 - 7) In (r V')] • (6) 
c 
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Fig. 3. The BFKL eigenvalue xil) ^-s a function ofj^ when 7i = 0, with 7 = 7i.+i7i. 
The single poles at 7 = 0, 1 correspond to the collinear (double logarithmic) limits. 
The saddle point for the Pomeron intercept occurs at = 1/2, while 7s = 0.372 
corresponds to the relevant anomalous dimension for saturation (see Sec. 6). 

In the above equation, ^ is a momentum scale related to the target, the 
integration contour C is parallel to the imaginary axis with < Re(7) < 1, 
and To (7) is proportional to the Mellin transformation of the initial cross 
section at y = 0. For instance, when we consider dipole-dipole scattering 
this initial condition is given by Eq. (B.7) in Appendix B and in this case 
ro(7) is equal to l/[2j'^{l-jf]. 

At high energies, i.e. when y — > 00, and with the dipole size r being 
fixed, one can perform a saddle point integration in Eq. (6). The sad- 
dle point of the x-function occurs at 71P = 1/2, with corresponding value 
x(l /2) = 4 In 2 as shown in Fig. 3, and then the dominant asymptotic energy 
dependence of the amplitude is given by 



where = 4asln2 is the so-called hard Pomeron intercept^. In Eq. (7), 

^ The expression for (T) in Eq. (7) vanishes in the strict large-A^c limit (qs = fixed, 
—i 0), because of its prefactor which is proportional to the initial value of the 
amplitude. But there is no real problem with that. First of all, and as we will 
see later on, the BFKL equation remains valid at finite-A^c- Furthermore, one can 
always assume a "modified" large-A'^c limit where all quantities, like the r.h.s. of 
Eq. (7), which are suppressed by factors, are taken as leading order effects 

provided they are enhanced by appropriate powers of the energy [6, 7]. For instance 
exp[2(x;py] is a leading effect, while exp[Li;py] is subdominant. 
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if represents the gluon density of the hadron, while n will be its dipole 
density if we assume that it is composed of dipolcs and this is the case 
in the large- A'c limit. The exponential increase is not surprising, since the 
BFKL equation is a linear one; in each step of the evolution, the change in 
the amplitude is proportional to its previous value. After all, this is what we 
had already expected from the summation of the series whose n-th term is 
given in Eq. (1). Notice also that, even though this is not explicitly written 
in Eq. (7) but easily inferred from Eq. (6), the dominant r-dcpcndcncc of 
the amplitude will be proportional to r. This is to be contrasted with the 
result of fixed order perturbation theory which is proportional to r^. The 
fact that the anomalous dimension differs from this perturbative result by 
a finite pure number is a unique characteristic feature of BFKL dynamics. 

The BFKL equation has also been solved in its most general case, 
i.e. without assuming rotational symmetry and for a given, but arbitrary, 
impact parameter [4]. In that case the system evolves "quasi-locally" in 
impact parameter space, and therefore the dominant energy dependence of 
the solution, as y ^ oo and with r and b fixed, is still given by Eq. (7), 
since the same eigenvalue dominates the corresponding integration. 

3. Pathologies of the BFKL Equation 

There are two major problems associated with the BFKL equation. The 
first is the violation of unitarity bounds. As we saw in Sec. 2 the amplitude 
at a fixed impact parameter increases exponentially with the rapidity Y , or 
equivalently as a positive power of the total energy s in the process, since 
Y = ln(s/.so) (the precise value of the scale sq is not important for our 
arguments). However, it should satisfy the unitarity bound^ 

r(r,6)<L (8) 

One should stress that this bound is different from the well-known Froissart 

bound [70-73], which states that any hadronic total cross section o"tot should 
not grow faster than the square of the logarithm of the energy, more precisely 

o-tot < ^ In^ s, (9) 

with the pion mass. Of course we should not expect to satisfy this bound 
by weak coupling methods. The BFKL equation, and all the equations that 
we will present later on in this article, will never treat properly the "edges" 

* This will be totally clear later on. See Eq. (24), where the precise definition of 
T(r, b) is given in terms of Wilson lines, thus including an arbitrary number of gluon 
exchanges with the target. 
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of the hadron, where long range forces become important and therefore the 
problem is genuinely non-perturbative. Nevertheless, there is no reason a 
priori which would imply that Eq. (8) cannot be fulfilled in perturbation 
theory. Notice that Eq. (8) is equivalent to the condition that the maximal 
allowed gluon density in QCD should be of order l/a^. Indeed, one has, 

^ a^a ~ ^ -77 ~ — , (10) 

where and a are gluonic creation and annihilation operators, while A 
is the gauge field associated with the wavefunction of the hadron. This 
maximal value of order l/a^ can be obtained in a heuristic way by setting 
the cubic and the quartic in A terms in the QCD Lagrangian to be of the 
same order. 

The second problem is the sensitivity to non-perturbative physics. As 

already said, the longitudinal momenta are strongly ordered in the course 
of the evolution. However the transverse coordinates of the dipoles (or 
equivalently the transverse momenta of the gluons) are not strongly ordered, 
and therefore the dipole kernel is non-local as clearly seen in Eq. (2)^. This 
non-locality results in a diffusion factor which accompanies the dominant 
asymptotic behavior given in Eq. (7), and which can be easily obtained by 
performing the Gaussian integration over 7 in Eq. (6) around the saddle 
point TP = 1/2. This factor reads 



with Dp = 2asx"(l/2) = 56asC{'i)- The form of the function suggests 
that it can be written as the solution to the one-dimensional diffusion equa- 
tion. That is, after the dominant exponential behavior has been isolated, the 
evolution can be viewed as a random walk in ln(r^/x^). Thus, independently 
of how small the initial dipole is, and after some critical value of rapidity, 
there will be diffusion to the infrared; dipoles with sizes bigger than Aq^j-, 
will be created and the weak coupling assumption will not be valid any 
more. One needs to emphasize that this is true even if we want to calculate 
the amplitude in the perturbative region. Since the diffusion equation is of 
second order, one can formulate the problem as a path integral from the 
initial dipole size to the final one. Then, clearly there are paths which go 

® In DGLAP [74] evolution, where one resums enhanced powers of ln((5^//i^) one 
encounters the "opposite" situation; the transverse momenta are strongly ordered, 
while the longitudinal ones are not. Therefore the corresponding DGLAP kernel is 
local in Q, but non-local in Y. 
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through the non-perturbative region. In this sense the BFKL equation is 

not self-consistent. 

Before moving to the next section, where we will discuss the solution to 
these problems, let us mention that the next to leading BFKL equation [27, 
28] shares the same pathological features. In that case one resums powers of 
the form as{asY)'^ , and this procedure gives rise to a contribution of order 
0{as) to the Pomeron intercept oop, as compared to the leading one. 

4. Unitarity, Saturation and Mergings of Pomerons 

Let us continue our heuristic approach and try to find out which is the 
basic element that is missing from the procedure we have followed so far. 
The first diagram in Fig. 4 shows one of the contributing diagrams to the 
BFKL equation. As indicated, this is of order asAYO{asip), where the 
factor as AY comes from the evolution step, the factor comes from the 
two couplings at the lowest part of the diagram and the factor ip, the target 
gluon density, is simply the upper part of the diagram. It is clear that the 
other two BFKL diagrams, which are not shown here, are of the same order. 
Following the same reasoning, the second diagram is of order asAYO{a1ip'^), 
since now there are four vertices at the lowest part, while one also probes the 
gluon pair density of the hadron. This diagram is suppressed with respect to 
the BFKL ones when (p <C l/cxs, i-e. at low densities, or equivalently when 
T <^ 1 and in that case it can be neglected. However it is equally important 
in the high density limit^, precisely in the region where the unitarity problem 
of the BFKL evolution appears. Therefore, following this "active" point of 
view in the projectile evolution as before, i.e. the r.h.s. of the equation 
contains what we have after the evolution step, we see that we should add 
to the BFKL equation the term which corresponds to the simultaneous 
scattering of the two dipoles {x,z) and {z,y), that is 



The resulting evolution equation is the first Balitsky equation [20]. There 
are a couple of points which need to be clarified in the above equation. 
The first is the notation "merge" for this particular term. This comes from 
the fact that the second diagram in Fig. 4 is equivalent to the third one. 
The latter corresponds to the evolution of the hadron in a "passive" point 
of view; now the red gluon is supposed to be emitted in the wavefunction 
of the target, and the two Pomerons, which existed before the evolution 

® In that case, this is a just typical diagram, since, in general, diagrams with more than 
two-gluon exchanges per dipole will be equally important. 
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Fig. 4. Estimation of diagrams contributing to the first Balitsky equation. 

step, merge to give rise to a contribution to the single dipole scattering 
amplitude. However, one should be very careful about the terminology and 
the interpretation of this merging process. When the hadron is boosted 
to higher and higher energies, there is no process which would lead to the 
"mechanical" recombination of the partons that it is composed of. The 
term merging corresponds to the 4^2 vertex connecting the upper part of 
the diagram, the hadron, and the lower part, the dipole. The second point 
that needs to be explained is the minus sign in Eq. (12). Intuitively, one 
would indeed expect this term to be negative, since, in almost any physical 
system, a merging process is supposed to tame the growth of the particle 
density. A better way to justify the sign, is to notice that we can express 
the first Balitsky equation in a more compact form, by introducing the 5- 
matrix Sxy = 1 — Txy (recall that T is the imaginary part of the scattering 
amplitude). Then we have 



which has a clear interpretation; the S-matrix for the dipole-pair to scatter 
off the target multiplied by the splitting probability gives the change in 
the S-matrix for the dipole-hadron scattering (after we subtract the term 
proportional to {Sxy) which gives the survival probability for the parent 
dipole). 

One can proceed and perform certain approximations to this equation. 
When the target hadron is a large nucleus, and for not too high energies, 
one can perform a sort of mean field approximation, by assuming that the 
two projectile dipoles scatter independently off the target, namely 



and similarly for the S'-matrices. Then one obtains a closed equation, the 
so-called BK (Balitsky, Kovchegov) equation [20, 21] (see also [75]). We 
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immediately notice that this equation has two fixed points; (i) (T) =0 4^ 
(S) = 1, which is an unstable fixed point, since no matter how small the 
initial amplitude is, it will start to grow and (ii) (T) = 1 '^=^ {S) = 0, which 
is the black-disk limit, and which is a stable fixed point for any perturbation 
in (T). Therefore the BK equation is much better behaved than the BFKL 
equation. It seems that the amplitude will never exceed the unitarity bound 
Tmax = 1, and the gluon density will saturate at a value of order 0{l/as). 
Furthermore the non-linear term cuts all the diffusion to the infrared [76, 77] 
and there is no sensitivity to non-perturbative physics any more. In fact, 
all diffusive paths that go beyond the saturation line (see below), will be 
eliminated by the non-linear evolution. 

Here it is appropriate to introduce the concept of the saturation momen- 
tum Qs- This is a line in the In{r'^fp) — Y plane (in general, the saturation 
momentum depends also on the impact parameter), along which the ampli- 
tude satisfies (T(r = l/Qs)) = const < 1. It is simply the border between 
the region where BFKL dynamics can be safely applied and the region where 
saturation has been reached and unitarity corrections have to be taken into 
account. At this moment, it is not hard to understand that more and 
more gluonic modes in the wavcfunction of the hadron will be saturated, 
as we keep increasing its energy. As we shall see later on in a detailed 
analysis, the saturation momentum increases exponentially with rapidity, 
i.e. ^ ^QCD exp(Asy). Thus, when the rapidity is large enough, one will 
have ct{Qs) 1, which means that the use of weak coupling techniques is 
justified. It is conceptually important to realize that this is a "QCD phase" 
where the physical system is dense but the coupling constant is weak. 

By no means we have given a rigorous derivation of the first Balitsky 
equation in this section. For example, when high density effects become sig- 
nificant, one should not restrict oneself to the two-gluon exchange approxi- 
mation (even though we didn't really make any use of it) when considering 
the scattering of a single dipolc with the hadron. Indeed, each coupling of 
either the quark or the antiquark of the dipole with the classical field A 
associated with the target is of order gA ~ y/Oisf, and thus in the high 
density limit (p ~ l/os all multi-gluon exchanges are equally important. 
Nevertheless, as we shall see in the next section, the first Balitsky equation 
remains as given here, even when we include multiple gluon exchanges. 

5. The Color Glass Condensate, the JIMWLK Equation and the 

Balitsky Hierarchy 

Perhaps the most elegant, modern and compete approach to describe 
the merging of Pomerons is the Color Glass Condensate (CGC) and its 
evolution according to a Renormalization Group Equation (RGE). This is 
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an effective theory within QCD, and its name is not accidental, since it 
corresponds to some of the basic features that it accommodates. Color 
stands for the color charge carried by the gluons, Glass stands for a clear 
separation of time scales between the "fast" and "slow" degrees of freedom 
in the wavcfunction, and Condensate stands for the high density of gluons 
which can reach values of order 0{l/as). 

The essential motivation for the formulation of the CGC is the separation 
of scales in the longitudinal momenta between the fast partons and the 
emitted soft gluons. Denoting by and /c^ the corresponding light-cone 
4- momenta, one has <C p^. This translates to an analogous separation 
in light cone energies k" = fc^/2fc+ p" = p^/2p+, which in turn leads 
to a separation of time-scales; the lifetime Ax^ ~ l/k~ of the soft gluons 
will be very short in comparison with the typical time scale ~ for the 
dynamics of the fast partons. Thus, even though the fast partons are virtual 
fluctuations in reality, they appear to the soft gluons as being a "source" 
which is x^-independent, i.e. as a static source. Furthermore, the source is 
random, since it corresponds to the color charge seen by the soft gluons at 
the short period of their virtual fluctuation; this happens at an arbitrary 
time and it is instantaneous compared to the lifetime of the source. The 
color charge density pa{x~,x) associated with this source at the scale p'^, 
propagates along the light cone x~ and the corresponding current has 
just a -|- component. This source is localized near the light cone within a 
small distance Ax~ ~ which is non-zero, but much smaller than the 

longitudinal extent ~ 1/k^ of the slow partons. One of the consequences 
is that the size of the hadron will extend in the longitudinal direction with 
increasing energy. 

Based on the above kinematical considerations, one can represent the 
fast color sources by a color current Ja{x~,x) = 8^^pa{x~,x), and the 
small-x gluons correspond to the color fields as determined by the Yang- 
Mills equation in the presence of this current, namely 



In principle, and in reality in certain gauges, one can solve this classical 
equation and obtain the gauge field A{p) as a function of a given source. 
Then, any observable O which is related to the field, for example the gluon 
occupation number or the amplitude for an external dipole to scatter off 
the hadron, can be expressed in terms of p, and its expectation value will 
be given by the functional integral 



{D,F-^^)^{x) = 8^^+Pa{x-,x). 



(15) 




(16) 
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It is obvious that Zy[p] serves as a weight measure and it gives the prob- 
abihty density to have a distribution p at a given rapidity Y (we assume 
that Zy[p] is normahzed to unity). This probabihstic interpretation for the 
source reUes on its randomness mentioned above; there is no quantum in- 
terference between different p configurations. In other words, one performs 
a classical calculation for a fixed configuration of the sources, and then 
one averages over all the possible configurations with a classical probability 
distribution. 

So far, in Eqs. (15) and (16) we have not used at all the small-,T QCD 
dynamics, and it is obvious that this must be encoded in the probability 
distribution Zy[p], which is the only quantity in Eq. (16) depending on 
rapidity. But before analyzing this aspect in more detail, we have to say 
that this idea to describe the soft modes of the hadronic wavefunction in 
terms of classical fields and probability densities was introduced in the MV 
(McLerran, Vcnugopalan) model [10], where a nucleus with large atomic 
number {A ^ 1) was considered. In this "static" model the only color 
sources are assumed to be the A x Nc valence quarks, which are taken to 
be uncorrelated for transverse separations such that |Aa;| < Aq^j^, so that 
the probability density is given by the Gaussian [10, 78] 



ZmvIp] ~ exp 



-1 



1 f'^-. p-{x)p- (x) 
2] ^"^I^ 



(17) 



where p? ~ K^^j^A^/^ is the average color charge density squared. Even 
though there is no y-depcndence in the model, a strong coherent color 
field of order 0{l/g) can be created due to the large number of nucleons 
and the solution of, the non-linear, Eq. (15) will lead to the saturation of 
the gluon density at a value of order 0{l/as), and with a saturation scale 
Ql{A)^Kl^^A^lHnA. 

Now let us see how and why the evolution of the probability distri- 
bution arises. Indeed, at rapidity Y , one should include in the source 
all those modes with longitudinal momenta (much) larger than A;"^, where 

Y = ln(l/x) = ln(P+/A;^") with P+ the total longitudinal momentum of 
the hadron. Now let us imagine that we increase the rapidity from Y to 

Y + Ay. Then, some modes that previously were "slow" , now they become 
"fast" and they need to be integrated over in order to be included in the 
source. This is a procedure which will crucially depend on the actual QCD 
dynamics and will lead to a change in the probability distribution. Still, we 
should emphasize that the theory at the new "scale" Y+l^Y is again defined 
through Eqs. (15) and (16), if we simply let F — > F -|- AF. A simple picto- 
rial interpretation of the evolution with Y of the probability distribution is 
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Y ^ ^ Y+AY ^ ^ Y+AY ' ^ Y+AY 

Fig. 5. Evolution of the probability distribution of color charges. 



shown in Fig. 5, where the red gluon at the rapidity Y' is a representative 
of the "semi-fast" modes which are integrated, i.e. Y <tiY' <^Y + AY. 

We shall not present here the derivation of this RGE, which is called the 
JIMWLK equation (even though a "sketch" of such a derivation is given 
in Appendix E). It can be found in the original papers [16-19] (see also 
[79,80]), while rather simple derivations can be found in [44,55,81]. The 
basic element is the resummation of Q;sln(l/x) enhanced contributions in 
the presence of a background color field. The latter is allowed to be strong in 
general, since one aims to describe properly the physics in the high density 
region. This equation can be given an elegant Hamiltonian formulation, and 
its most compact representation arises when it is expressed in terms of the 
color field a{x^,x) = A'^{x~,x) in the Coulomb gauge. It reads 



H 



a, — 
oa 



ZY[a], 



where the explicit form of the Hamiltonian is [17, 18, 26] 



H 



167r3 



M-, 



ab § 



(18) 



(19) 



The dipole kernel^ is readily recognized, while the Wilson lines V in the 
adjoint representation arise from the gluon propagator of the integrated 
modes and they are given by 



a 



P exp 



ig / dx-a''{x-,x)T'' 



(20) 



The longitudinal coordinate in the functional derivatives is to be taken at 
oo, since x~ ~ l/A:"*" ~ e^/P+ oo; as we pointed out earlier, the hadron 



^ The dipole kernel appears in the above equation when one assumes that the Hamil- 
tonian will act on the observables defined below in Eq. (21). In general one should 
replace the dipole kernel Mu-oz by 2{u — z)-{z — v)/{u — z)'^{z — v)^. 
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extends in the longitudinal direction, and all the "action" is expected to take 
place at the last layer of rapidity. Finally, the two functional derivatives 
appearing in Eq. (19) correspond precisely to the two outgoing legs in the 
bottom of the diagrams in Fig. 5. The appearance of only two such legs is 
related to the certain class of diagrams which are effectively resummed by 
the JIMWLK equation; during a single evolution step, only the two point 
correlation function (pp) changes, and this comes from the fact that the 
color field, or equivalently the sources, are assumed to have values (much) 
larger than g. This is certainly a well-defined approximation, but with a 
decisive influence on the outcomes of the theory. We shall not discuss more 
details here, but we shall return to this issue and to its proper treatment, 
in the next sections. 

Given a Hamiltonian, the natural question that one may ask is which 
the observables are. They can be nothing else than the quantities already 
appearing in the Hamiltonian, and to be more precise, they will be gauge 
invariant operators built from Wilson lines. For our purposes, we shall only 
consider Wilson lines in the fundamental representation and then the most 
generic form of such an operator will read 



Notice that a Wilson line has a direct physical interpretation [82]. Let us 
consider the scattering of a left-moving quark off the classical field created 
by the (right-moving) hadron. The general expression for the S'-matrix in 
the interaction picture is 



where T stands for a time ordered product, and the relevant part of the 
QCD Hamiltonian for this process is Hi = —g^^~ A^Tp. Here we have used 
the fact that the quark is a fast left mover and therefore the 7~ matrix 
dominates the inner product. Now, since the trajectory is "eikonal", the 
transverse and -|- coordinates of the particle are fixed, and they are chosen 
to be x' = X, and (by convention) x'~^ = x+ = 0. Therefore, so long as the 
color independent part is concerned, one can make the substitution 



Since X is increasing along the trajectory, one can replace the T-product by 
a Path ordered product. Therefore the ^-matrix is equal to the Wilson line 
Vx given in Eq. (20), but in the fundamental representation, i.e. one needs 



0[a] = tT{vlVx,VlVx, . . ■)tr{VlVy, ...)... 



(21) 




(22) 



^{x') j-A+{x') il^{x') 5{x' - x) 5{x'+ )A+{x). 



(23) 
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to do the replacement t"". Notice that all kind of multiple exchanges 

are included in this procedure, and they can be recovered by expanding the 
Wilson lines to a given order in the coupling constant. It is not hard to 
understand that the ^-matrix for a dipole [x, y) to scatter off the target 
will be given by the gauge invariant expression 



Sxy = — tic{V^Vy) = l-Txy = l-{-{al- alf + 0{g\ (24) 



where we have also expanded to second order to obtain the scattering am- 
plitude in the two-gluon exchange approximation for later convenience, and 
with the field a" in this expansion being integrated over x~ . 

As discussed earlier, the evolution of the expectation value of a generic 
operator O will come from the evolution of the probability distribution 
Zy[p\. Using Eqs. (16), (18) and (19), and after a functional integration by 
parts we find 



Now we are ready to write the equations obeyed by the scattering ampli- 
tudes. Naturally, we start from the scattering of a single projectile dipole. 
Using Eqs. (19), (24) and (25), and as we show explicitly in Appendix F, we 
arrive at the first Balitsky equation which we rewrite here for convenience 



Notice that the above equation is valid for a finite number of colors and the 
same will be true for the BFKL equation which arises in the limit (T) = 
1— (S) <C 1. As we have already discussed, this is not a closed equation and 
one needs to find how {SxzSzy) evolves with rapidity. Following the same 
procedure, and as we show again in Appendix F, one obtains the second 
Balitsky equation, which reads 




(25) 




(26) 



z 




w 




w 




w 
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where we have defined the "quadrupole" operator 

1 



Qxzwy 



No 



^^V^y^V^zVyVlV,) -ii{V^Vy) . (28) 



At this point a few comments need to fohow. We immediately notice that 
the first two terms in the second Bahtsky equation can be obtained just 
by applying the Leibnitz differentiation rule to the first Balitsky equation. 
This has a natural interpretation in terms of projectile evolution as the 
two dipoles (a;, z) and {z, y) can evolve independently; cither of these two 
dipoles can split into two new dipoles which subsequently scatter off the 
unevolved target hadron. However, there is the third term which goes be- 
yond this simple dipolar evolution. This is not surprising since the small-x 
gluon emitted by the one of the two original dipoles can be absorbed by 
the other dipole (more precisely, emitted by the other dipole in the complex 
conjugate amplitude when calculating the wavefunction of the projectile) 
and therefore no dipole survives after the evolution step. Thus, the pro- 
jectile system is lead to a more complicated multipolar state, which can be 
naturally called a color quadrupole because of its dependence on four trans- 
verse coordinates. For the next step one will have to write not only the 
evolution equation for the thrcc-dipolcs operator, but also the one for this 
quadrupole operator. It becomes clear that the complexity in the structure 
of this hierarchy of equations, which is called the Balitsky hierarchy, will 
rapidly increase as we proceed to describe the evolution of "higher-point" 
functions. Following this "active" point of view in the projectile evolution, 
we see that is wavefunction will be successively composed of 

one dipole two dipoles three dipoles + quadrupole 
— > . . . ^ dipoles -|- higher multipolar states. 

Furthermore, we should emphasize that no factorization like the one in 
Eq. (14) will solve the infinity hierarchy. Nevertheless, we see that the 
two fixed points of the BK equation, when generalized properly, satisfy this 
system of equations, (i) When the classical color field a" vanishes there is 
no scattering; the expectation values of all the Wilson lines become equal to 
unity, i.e. {S) = {SS) = {Q) = ■ ■ ■ = 1 and (ii) when the color field becomes 
large, we reach the black-disk limit; the Wilson lines oscillate rapidly [83] 
and their expectation values vanish; i.e. {S) = (SS) = (Q) = • • • = 0. 

One expects drastic simplifications in the large- iVc limit where the de- 
grees of freedom can be chosen to be the color dipoles. Indeed, the last 
(quadrupole) term in Eq. (27) is of order 0{1/N^), negligible at large-A^c 
when compared to the first two (dipolar) terms which are of order 0{1) and 
therefore the evolution will proceed only through dipolar states. Now one 
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can see that the hierarchy becomes consistent with a factorization of the 
form [84, 85] 

{Sxiyi ■ ■ ■ Sx^y^) = C {Sxiyi) ■ ■ ■ {^XkUk) ^ (^9) 

with c an arbitrary constant. Then the whole large- AT^ hierarchy collapses 
to a single equation, the BK equation but with a coefficient c in front of the 
{Sxz){Szy) term. The asymptotic fixed point of this equation is clearly 1/c, 
but it is not clear whether a value c 7^ 1 has a simple physical interpretation 
or not. In any case, if the initial conditions at a rapidity Yq arc of factorized 
form, this factorization will be preserved by the large-iVc evolution. But 
even if the initial conditions are not of this type, no new "correlations" 
will be generated by this large-iVc evolution; only the initial ones will be 
propagated to higher rapidities. 



6. The Saturation Momentum 

With the theory of the Color Glass Condensate being theoretically es- 
tablished, one of the central problems has been the determination of the 
saturation momentum Qg which, as we discussed at the end of Sec. 4, is the 
momentum scale at which we start to approach unitarity limits. In principle 
one should be able to calculate the energy (rapidity) dependence of Qs (at 
least asymptotically), since all the dynamics is contained in the JIMWLK 
equation, however its precise value cannot be determined since it will depend 
on initial conditions and thus on details of non-perturbative physics. Recall 
that the saturation momentum is an intrinsic property of the target hadron, 
since it is also the scale where the gluon density saturates. And the only 
scale associated with the hadron is Aqcd- Given the above considerations, 
and the fact that the pure BFKL evolution leads to an exponential growth, 
one may guess that the saturation momentum will be of the form (for fixed 
coupling, neglecting the impact parameter dependence and at asymptotic 
energies) 

Q's^cAla^"-^^, (30) 

where the constants A and (3 should be calculable, while the constant c would 
require the knowledge of the non-perturbative structure of the hadron. 

Even though the problem seems difficult at a first sight due to the com- 
plicated structure of the Balistky-JIMWLK equations, one can do certain 
logical simplifications that render the calculation tractable. As we saw, the 
hierarchy reduces to the much simpler BK equation in the large- A^c limit 
and under a mean-field approximation, and then, presumably, the most that 
one may lose is corrections. Furthermore, one can in fact use BFKL 
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Fig. 6. Scattering amplitude and saturation momentum in the logarithmic plane. 

dynamics, if one is careful enough to treat properly the effects of the linear 
terms. In turn, this means that we may not even lose the possible 
corrections, since the BFKL equation is valid at finite- iVc. 

Now let us have a look at Fig. 6 where the expected qualitative behavior 
of the scattering amplitude in the logarithmic plane ln(Q^/^^) — y is shown. 
Here Q should be thought as the inverse of the dipole size or as the transverse 
momentum of a gluon in the hadron. Along Line-1, which corresponds to an 
anomalous dimension 7 — > 0, the amplitude decreases since the momentum 
increases too fast. This corresponds to DGLAP evolution (in the double 
logarithmic limit) where one resums powers of as ln(Q^///^) (x ln(l/x)). In 
this case there is again a cascade of partons inside the hadron, but these 
partons are very small in size and they will never overlap to form a high 
density system. Along Line-2, which corresponds to an anomalous dimen- 
sion 7p = 1/2, the amplitude increases, since the energy increases while the 
momentum remains fixed. This is the hard Pomeron intercept line and cor- 
responds to the saddle point of the BFKL eigenvalue function. Starting from 
a value of order O(a^) and after covering a rapidity interval ~ ln(l/a^)/a;p. 
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the amplitude will eventually "hit" the unitarity /saturation region. 

Therefore, there must be some "critical" line between those two lines 
which belongs to the region of linear evolution and along which the ampli- 
tude remains constant, for example of order 0{as) <C 1. Clearly this line 
must correspond to an anomalous dimension, say 7^, such that < 7^ < 1/2. 
Then the saturation line which corresponds to constant amplitude of order 
0(1) but smaller than 1, for example 1/2, will be parallel to the critical line, 
and therefore characterized by the same anomalous dimension and energy 
dependence. The reason why the two lines are parallel, may not be so clear 
at the moment, but we will try to justify it in a while. 

The question is whether or not we can use the BFKL dynamics to de- 
termine this critical line. Even though the line belongs to the linear region, 
the answer is negative if one wants to get the prefactors correct, i.e. the 
value of P in Eq. (30). The reason is that as the system evolves along this 
line, and after we have isolated the leading exponential behavior, there we 
will be some paths, really in the functional integral sense, that go through 
the saturation region and then return to the linear one. These diffusive 
paths are absent when the full BK equation is equation is considered, and 
therefore should be cut. This can be done by using an absorptive boundary 
just beyond the saturation line, which will mimic the effects of the non- 
linear terms for the problem under consideration. Put it another way, this 
procedure is equivalent to a self-consistent solution of the BFKL equation 
with the boundary condition {T{Q = Qs{Y))) = c < 1. 

Now we are ready to determine the "slope" of the critical, and therefore 
of the saturation, line. We impose two conditions for the exponent in the 
solution to the BFKL equation as given in Eq. (6)^. First we require the 
exponent to be zero so that the amplitude be constant. Second we impose 
a saddle point condition, which is a valid approximation in the asymptotic 
limit agY ^ 1. Then we find that the anomalous dimension is determined 
by the saddle-point of x{l)/{^~l) is given by [2, 77, 86-88] 

X(7.) + (1-7.) X' (7.) =0^7. = 0.372. (31) 

The energy dependence of the saturation momentum is better expressed 
in terms of its logarithmic derivative which reads [77, 88] (with the leading 
term already known from [86, 87]) 

_ dlnQ^ _ x(7.) 3 1 2.39 



For our convenience, here we consider a case where the ampUtude is averaged over 
impact parameter and therefore it is dimensionless. Then it is still given by Eq. (6) 
but the factor 2tv/ij? is not present. 
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while for the scattering amplitude one obtains [77, 88] (and with the absence 
of the logarithmic modification, already known from [87]) 



DsY 



(33) 



a form which is valid in the region <C <C Q^g ex.\){D sY) , and where 
the diffusion coefficient is Dg = 2asx"{'ls) = Q^ctg. 

Now let us discuss the results. As we see in Eq. (31), whose graphical 
solution is shown in Fig. 3, the relevant value of the anomalous dimension 
for saturation lies indeed in the interval (0, 1/2). We should say that the 
eigenvalue xils) will be selected, so long as the initial condition contains the 
corresponding eigenfunction and this will be true for all interesting cases. 
Notice that 7^ is a pure number, and this would have never been obtained 
by applying pure DGLAP evolution. The latter always gives anomalous 
dimensions which start at order 0{as). 

In Eq. (32) we see that the leading contribution to the "intercept" of 
the saturation line is totally fixed by BFKL dynamics. It is not too difficult 
to see how the subleading correction arises, and to this end let us go a few 
steps back in the derivation of Eqs. (32) and (33). After we have performed 
the Gaussian integration around the saddle point 75, the solution reads 
T oc {Qo/Q^)^~'^^'4's, with Qo containing only the leading behavior of Qg, 
and where i/jg satisfies the diffusion equation. The "standard" solution to the 
diffusion equation behaves as Y^/"^ (times the exponential diffusion factor), 
but in the presence of an absorptive boundary the survival probability of 
the "particle" becomes smaller and tpg is proportional to l/F^/^. Thus, 
by combining this prefactor and the leading behavior Qq, one obtains the 
correction written in Eq. (32). Notice that the saturation effects lead to a 
slower increase of the saturation momentum, as they should. Even though 
the second term vanishes when y ^ 00, it cannot be neglected, since upon 
integration of Eq. (32) it will generate a y-dependent prefactor in Qg. This 
1/Y term in the logarithmic derivative of the saturation momentum should 
be interpreted as l/R'jy, where Rd ~ VY is the diffusion radius. This radius 
is practically the available phase space in logarithmic units of transverse 
momentum. In reality this phase space extends to infinity, since there is no 
boundary to the ultraviolet, but in practice it is only the space inside the 
diffusion radius which will contribute to the subsequent steps of evolution. 

Finally, as we see in Eq. (33), the scattering amplitude has a scaling 
form [77,87], so long as lio.{Q'^/Qj) <C VDgY so that the exponential dif- 
fusion factor can be set equal to unity. That is, the amplitude does not 
depend on and Y separately, but only through the variable Q^/Q^- The 
pure power is not an unexpected result; the BFKL evolution generates the 
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anomalous dimension 7^ which modifies the behavior ~ of fixed order 
perturbation theory. The logarithm is generated by the absorptive bound- 
ary when solving the diffusion equation. Notice that both terms in Eq. (33), 
the power and the power modified by the logarithm, are exact degenerate 
solutions to the BFKL equation [89], and thus the final solution is just their 
linear combination. Such a scaling behavior, which will be preserved even in 
the running coupling case (but in a more narrow window in Q^), but which 
will be violated by the evolution equations that we will discuss later on in 
Sec. 10, is consistent with fits [57, 58, 60] of the low- a; data in DIS. 

Here it is appropriate to mention that the solution of the BK equation 
close to the unitarity limit Aq^j-, <C ^ is also of scaling form and as 
we show in Appendix G it reads [37, 83, 90] 



(5) = 1 - (T) « exp 



7. ^Q: 



2x(7.) 



(34) 



Thus, one naturally expects the scattering amplitude to satisfy scaling ev- 
erywhere from deep inside the saturation region up to momenta which be- 
long in the region of linear evolution. This explains why the critical and 
saturation lines are parallel to each other. 

Eqs. (31), (32) and (33) have been confirmed by studying the anal- 
ogy to the (nonlinear) FKPP (Fisher, Kolmogorov, Petrovsky, Piscounov) 
equation [88,91]. Furthermore, the last asymptotic term of As(F) which is 
independent of the initial conditions, and whose behavior is ~ l/Y'^/^ has 
been obtained in the same fashion in [92] . The full JIMWLK hierarchy has 
been solved numerically on the lattice [93] and the results agree with the 
ones we presented in this section. Moreover it was found that violations of 
factorization in the scattering amplitude correlations are extremely small. 

Before closing this section, let us see how these results are modified 
when we consider BFKL dynamics at the next to leading level. The calcu- 
lation of the next to leading order (NLO) correction to the BFKL kernel 
was completed in [27,28]. However, this negative correction turned out to 
be larger in magnitude than the leading contribution for reasonable values 
of the coupling, say ois ~ 0.25. Even worse, when cxg > 0.05 the full ker- 
nel has two complex saddle points which lead to oscillatory cross sections 
[94]. But it was immediately recognized that these large corrections emerge 
from the collinearly enhanced physical contributions [95-98] . A method was 
developed to rcsum collinear effects to all orders in a systematic way and 
the resulting renormalization group (RG) improved BFKL equation was 
consistent with the leading order DGLAP [74] equation by construction. 

But before going to the NLO case, let us consider the leading BFKL 
kernel, but with a running coupling ^^(Q^). One expects two major changes 
with respect to the fixed coupling case. Since the saturation momentum 
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increases with rapidity, the coupling will decrease as we evolve close to the 
saturation line, and Qg will increase much slower. Moreover, the integration 
of the quadratic fluctuations around Qg will be affected by the fact that the 
coupling is not a constant quantity, and therefore the mechanism of diffusion 
will be modified. An analytic expression can be given and it reads [77, 89] 
(with the leading terms already known from [2, 87]) 

Xs = J^('--^^4j,)=^^-. '-''l. (35) 



(36) 



6(1-7,) Vr 4 r5/3; Vy + Yo {Y + Yo)^/^ 
for the logarithmic derivative of the saturation momentum, and 

,2^l-^s , ,..^2/^2^ , .X r 21n^(Q2/g2) 

3x"(7.)r J 



for the scattering amplitude. Here we have defined b = {llNc—2Nf) j [YINq)^ 

T = ^2x{7s){Y + Yo)/[b{l-ls)], Dr = {x"{ls)/[2x{ls)]V/' = 1-99, Ai is 
the Airy function and = —2.33 is the location of its leftmost zero. The 
second expression for A, corresponds to ATj = 3 flavors. Notice that the 
flrst term in A, is equal to the flxed coupling result, in the sense that it 
may be written as x{'ys)(is{Q'i) / {^~7s) ■ The second term is negative since 
it accounts for the contribution of the boundary (and the prefactors) and 

it has a parametric form a^^^ [89], a well-known type of correction in NLO 
BFKL dynamics [97,99]. We should notice that in this case, Qg will be 
always proportional to Aqcd at high rapidities 1" » Yq [77, 89, 100]. This is 
in contrast to the fixed coupling case where Qg is proportional to the initial 
saturation scale, if such one exists. For example if the target hadron is a 
large nucleus, the square of the initial saturation momentum and thus the 
constant c in Eq. (30) is proportional to A^^^ In A. Finally we note that the 
expression for the scattering amplitude takes a scaling form in the window 
lii{Q'^/Ql) <^ Dj.T^/^, and this form is exactly the one found in the fixed 
coupling case. 

Even though one cannot have such nice analytic expressions when con- 
sidering the improved kernels, the results are not so hard to obtain since 
they involve the numerical solution of algebraic transcendental equations 
[89]. We summarize all cases in the plots in Fig. 7. In order to understand 
properly all the results, we have also plotted the first three lines which cor- 
respond to the analytic expression given in Eq. (35). Line-a is the first term 
with lo = while Line-b corresponds to the same term but with a typical 
value for Yq which is of order 0(1). Line-c stands for the full expression in 
(35). Line-d and Line-e represent the improved kernels at leading and next 
to leading order respectively. Notice that all the lines will "merge" together 
at very high values of rapidity, since the leading term will become more 
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Fig. 7. The logarithmic derivative of the saturation momentum As as a [unction of 
the rapidity Y for running coupling. 



and more dominant over the corrections as the couphng decreases along 
Qg. The NLO result is stable in the sense that it adds a small correction 
to the leading one. We should mention that the NLO result is practically 
constant, ~ 0.3 in the region y = 5 — 10. This is in good agreement 
with the phenomenology [57-60]. However, it is not clear whether this the- 
oretical value should be trusted or not. As we shall see in the forthcoming 
sections, the JIMWLK equation misses important contributions, which re- 
sult in corrections that are much more significant than any NLO BFKL 
correction. 



7. Deficiencies of the Balitsky-JIMWLK Hierarchy 

Even though the Balitsky-JIMWLK hierarchy encompasses nicely the 
BFKL evolution and the merging of Pomerons at finite- A'^c, it faces certain 
crucial problems. These problems, which we will immediately discuss, are 
related to each other, and therefore it seems that some unique element is 
missing from the CGC effective theory in the form it has been developed so 
far. 

(i) The first problem is the extreme sensitivity of the JIMWLK equation 
to the ultraviolet. Since in the high momentum region the non-linearities 
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are unimportant, wc understand that this issue can be analyzed within the 
BFKL evolution. Let us try to reconstruct the solution in two (or more) 
global steps by completeness. To be more precise, let us first assume that 
we find the solution T{Q,Y) by evolving the system from zero rapidity to 
rapidity Y. Now let us imagine that we evolve from zero to, say, Y/2. 
Then we can consider the solution T{Q, Y/2) as the initial distribution and 
subsequently evolve to Y to obtain T{Q, Y). As we show in Appendix H, the 
solution obtained from this procedure will coincide with the one obtained 
from the single global evolution step, so long as we include (at least) the 
contributions from all momenta such that ln{Q'^ /Q^) < ^JD^YJ^, in the 
initial condition at y/2. Of course there is no reason to cut the momenta 
that lie outside the diffusion radius, but this algorithm reveals the width of 
the momentum phase space which is important for a self-consistent solution. 
This feature brings us in a quite embarrassing situation; as Y increases, this 
phase space will open up more and more to momenta above the saturation 
line, and moreover the big numerical value of the coefficient Ds will make the 
problem even worse. For instance, when one finds the saturation momentum 
to be a few GeV, at the same time one is sensitive to momenta a few orders of 
magnitudes above. Notice, that this explains why in the numerical solutions 
of both the BK [64, 76, 101-104] and the JIMWLK equation [93], one had 
to go very far to the ultraviolet in order to obtain a reasonably accurate 
solution. In the running coupling case the situation is somewhat better, 
since the coupling decreases at higher momenta, and thus the effects of 
these seemingly non-physical contributions are reduced. Indeed, as we saw 
in Sec. 6, the "diffusion" radius increases much slower, namely Rb ~ Y^l^ . 
Nevertheless, the theoretical problem still exists. 

(ii) Quite surprisingly, the second problem is the violation of unitarity. 
Here we shall not analyze the argument in detail [34], but only indicate 
its essence. Assume that we want to calculate the amplitude close to, but 
above, the saturation line in the two ways we described in the previous 
paragraph. Then one will have 

i>c = r~^rir2, (s?) 

where T\ and Ti denote the contributions of the two successive steps. It 
is clear that for T\ < the above equation imposes that the second step 
satisfy T2 > 1. Thus, all the paths which go through the region to the 
right of the fluctuation line (the terminology will be understood in a while) 
in Fig. 6 will violate unitarity in the intermediate steps of the evolution. 
Returning to the problem we discussed in (i) , and noticing that the diffusion 
radius extends to the region where the amplitude can be much smaller than 
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a^, wc see that these contributions from the ultraviolet region must be 
indeed non-physical. 

These unitarity violating paths were "discovered" when dipole-dipole 
scattering was studied in the presence of saturation with the additional, 
but clearly natural, condition that the amplitude be Lorentz invariant [34], 
even though at that time it was not realized that they are part of JIMWLK 
evolution. A calculation of the saturation momentum was performed by 
cutting these paths with an absorptive UV boundary (in addition to the IR 
one), and the corrections found were huge. In Sec. 15 we shall review/obtain 
this result in a very simple way. We should say here, that in [35] a spec- 
tacular analogy with the physics and the results of the stochastic FKPP 
(sFKPP) equation was observed, and the significance of fluctuation effects 
in the low density region T ~ '^=> </? ~ was understood. Still, the fact 
that one needs to go beyond the JIMWLK equation was not yet realized. 

(iii) The third problem in the evolution equations we have discussed 
so far, is the absence of Pomeron splittings [38]. As we have seen, the 
mechanism of Pomeron mergings, which was essential to describe properly 
the physics near the unitarity limit, is described by the JIMWLK equation. 
If we expand the Hamiltonian in powers of g, then each term will involve at 
least two factors of the color field a and exactly two functional derivatives 
with respect to a. Then a typical term in the evolution equation of the 
n-th point correlator of the color fields will have the structure (suppressing 
the coordinates' dependence, which is not important for the forthcoming 
argument) 



So, as we already knew, the JIMWLK Hamiltonian can describe BFKL 

dynamics and Pomeron mergings. But the natural question at this point 
is, "how could we have two or more ladders in the first place?" Clearly 
JIMWLK cannot do that, since one would need m < ra in Eq. (38). One 
option would be to consider a large nucleus target, which contains many 
valence quarks and antiquarks. These sources can evolve with rapidity and 
produce many BFKL Pomerons, which will merge when saturation becomes 
important. However, this is just a special case because of its particular ini- 
tial condition, and the dynamical problem is not solved. Furthermore, even 
in the nucleus, there will always be some dilute "tail" corresponding to the 
high momentum modes. After some evolution, and since the saturation 
momentum increases, these modes will need to saturate. But still there is 
no dynamics to produce the corresponding Pomerons which will eventually 
merge. Thus, the only solution to this problem is to find how QCD will 




m>n. (38) 
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give rise to the Pomeron splittings. Tlicn indeed, one could start, for ex- 
ample, even from a single bare dipole, and end up with a fully saturated 
wavefunction. As we shall see in the following sections, when "completing" 
the theory by including the diagrams which were "forgotten" [38,39], and 
which lead to the splittings of Pomerons, we will also automatically solve 
the two problems presented in (i) and (ii). 

8. The Missing Diagrams 

Since one of the basic problems of the JIMWLK equation is the ab- 
sence of Pomeron Splittings, one should consider diagrams like the third 
one in Fig. 8 in order to resolve the issue. Indeed, this diagram corre- 
sponds to a transition from one to two Pomerons. If one wants to "measure" 
two Pomerons, and in the two-gluon exchange approximation, one needs to 
probe the hadron wave function with two projectile dipoles, and this is what 
is shown in the typical diagrams in Fig. 8. The first diagram corresponds to 
normal BFKL evolution, with the one of the Pomerons being a spectator. 
As indicated, it is of order asAYO{a'^ip'^); a^Ay for the evolution step, 
for the four fermion-gluon vertices and cp'^ since there were two Pomerons 
before the step. Similarly the second diagram, which corresponds to a 2 ^ 1 
Pomeron transition, with the third Pomeron being a spectator, is of order 
asAYO{a^ip^). Both diagrams are described by the JIMWLK equation^. 
The third diagram is of order a^AYOia^ip), since there are six vertices 
and one Pomeron before the step. Clearly this diagram is suppressed with 
respect to at least one of the first two diagrams when » a^, and this is 
the reason why it was neglected in the derivation of the JIMWLK equa- 
tion, which aimed to describe a high density system. However, due to the 
non-locality of the evolution kernel which leads to the ultraviolet diffusion, 
this diagram will give significant contributions through the intermediate 
evolution steps. It is rather crucial to notice that this diagram becomes 
important when (yj ~ ofs or equivalently when T ~ a^, which corresponds 
precisely to the line at which JIMWLK faces its unitarity problem. This is 
the fluctuation line we have already shown in Fig. 6, and which divides the 
region of normal BFKL evolution from the low density region. The latter 
is characterized by fluctuations and this explains in a way why we need (at 
least) two projectiles dipoles. By probing with one dipole we will measure 
only the "one-point" function, say the gluon occupation number. But this 
will not release any information about the two-point function, since in a di- 
lute system the pair density is affected by the fluctuations and is not simply 

® Also diagrams like the first one, but with the soft gluon connecting the two ladders are 
included; these diagrams are of the same order in Us and (p, but they are suppressed 
at large Nc. 
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a,AYO{al^^) a,AYO{al^^) a,AYO{al^) 

Fig. 8. Diagrams contributing to the evolution of two BFKL Pomerons. The third 
diagram was "missed" in JIMWLK. 



the square of the single density. Thus one needs to probe with two dipoles 
in order to measure the non-trivial low density correlations. This analysis 
also implies that, presumably, the first Balitsky equation will not change. 
But of course, any mean field approximation to this equation will not be 
valid any more, since the presence of fluctuations will have a significant im- 
pact on its non-linear term. In the next section we shall shortly describe 
some elementary but essential, for our purposes, features of the color dipole 
picture. We shall return to give a more quantitative analysis of what we 
discussed here, in Sec. 10. 



9. Evolution of Dipoles 

Since we have realized that what we need to correct in the JIMWLK 
equation is its low-density limit, let us look more carefully at the non- 
saturated part of the wavefunction of a hadron. We will first consider the 
large- A'^c limit case, that is, the hadron is supposed to be composed of color 
dipoles. When the dipole density is not too high, i.e when n <C l/ttg, the 
emission of a soft gluon from a dipole will not be affected by the remain- 
ing surrounding dipoles. We have already written the (positive and well- 
defined) probability for this emission in Eq. (2) with the derivation given in 
Appendix A. From now on we shall not present the original formulation of 
this dipole picture [5-8], rather we will follow the procedure developed in 
[105] which is better suited to our purposes. 

One can write a master equation describing the evolution of the probabil- 
ities to find a given configuration. To be more specific, a given configuration 
is characterized by the number of dipoles and by — 1 transverse coordi- 
nates {z^} = {zi, Z2, . . . , zj\f-i}, such that the coordinates of the N dipoles 
are {zq,zi), {zi, Z2),. ■ ■ ,{zn-i, zn), with zq = uq and zn = vq, assuming 
that the initial state was a dipole {uq,vq). The probability P/v({zj};y) to 
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find a given configuration at rapidity Y obeys 

oTT 2^M[Zi-i, Zi+i, Zi)}^N-l{Zi, ZN-l,y ) 



dY 27r 

1=1 

TV 

a 



^ I M{z^-uz^,z)Pf^{zu ...,zi^-v^), (39) 
1=1 , 



where a slashed variable is omitted. The interpretation is quite obvious; 
while the gain (first) term describes the formation of an A^-dipoles state 
through the splitting of a dipole in a pre-existing state with A'^— 1 dipolcs, the 
second term describes the emission of a soft gluon from the A''-dipoles state, 
which leads to a state with A^^+l dipoles and therefore to a loss of probability. 
It is straightforward to show that the total probability is conserved, which 
is one of the requirements that Eq. (39) be well-defined. The expectation 
value of an operator O which depends only on dipole coordinates is given 
by 

oo „ 

(0(y)) = E / ^^nPn{Xz,}-Y)0^{{z,}\ (40) 

where the phase space integration is simply dF^v = d^Zid^Z2 . . . d^Ziv-i- 
Then by using the master equation (39) one can show that 

a(o(y)) _ a °° " ^ 



dY 

N=l 



^Yl dTr,PNi{z^}■,Y)Y I M{z,^^,z,,z) 

x[-0^({zj) + 0jv+i({^.,^})], (41) 



where the z argument in Ca +i is to be placed between and Zj. In what 
follows, we shall use Eq. (41) to derive evolution equations for the dipole 
number densities. Consider first the average number density of dipoles at 
(«, v). The corresponding part for an iV-dipole configuration is 

N 

nAr(M,'y) = E^(zj_i - u)5(zj - u), (42) 

.7=1 

so that {uuT^ will be given by the r.h.s. of Eq. (40) if we let Oat tin. 
Then by making use of Eq. (41) and after relatively simple manipulations 
one arrives at the evolution equation for the average of the dipole number 
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dY 



Fig. 9. Evolution of the dipole density for a dilute hadron at large Nc 

density {riuv) which reads 
d{nuv) a. 



dY 27r 

= ^ 



uzv )-M uvz {f^uv)] 



(43) 



This is simply the BFKL equation for the dipole density, and its pictorial 
representation is given in Fig. 9. This equation should be read in the "pas- 
sive" point of view, in contrast to Eq. (3). That is, as it was obvious in 
the derivation, the r.h.s. contains what was existing before the evolution 
step. The first term is proportional to the probability for a dipole (w, z) 
to split into two new dipoles {u^v) and {v,z) times the initial density at 
[u,z). Since we want to find the change in the density at {u,v) only the 
first of these two child dipoles is measured. Similarly for the second term. 
The last term is proportional to the probability for a dipole (w, v) to split 
into two new dipoles times the initial density at (it, v) and naturally gives 
a negative contribution. Notice that the loss term in Eq. (39) was crucial in 
order to obtain this virtual contribution. As we show in Appendix C, this 
equation leads to the BFKL equation, Eq. (3), for the scattering amplitude 
of a projectile dipole off the dilute hadron, when we assume the two gluon 
exchange approximation which is appropriate in this dilute limit. 

Now we can follow the same procedure to derive the analogous equation 
for the dipole pair density [38]. The corresponding part for a given A^-dipole 
configuration is 



,(2)/ 



N 



{uiVi;u2,V2) 



-ui)5{zj -vi)5{zk-i-U2)5{zk-V2), (44) 



j,k=l 
37^k 



and following the same procedure as in the case of the dipole density, we 
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dY 




+ 2 



Fig. 10. Evolution of the dipole pair density for a dilute hadron at large Nc- The 
loss term of the normal BFKL evolution is not shown. 



arrive at the evolution equation 



dY 



as 
2tt 



UlVlZ 



(8) {n. 



(2) 



uivi;U2V2 / 



+ 1^2. 



(45) 



with the notation introduced in Eq. (43). The two terms on the r.h.s. cor- 
respond to the two typical contributing diagrams in Fig. 10. The first term 
(which in turn is a sum of three terms) corresponds to the BFKL evolution 
of the first dipole while the second dipole remains a spectator. The second 
term corresponds to a single dipole initial density and the mother dipole 
splits into two new dipoles both of which are "measured" . We shall refer to 
this term as the "splitting" term, and this in the sense that a lower moment 
of the density gives rise to a higher moment of the density. 

Of course one can continue and write the evolution equation for the 
K-th density [106]. We shall not do it here, since it is just a matter of 
proper combinatorics. It is clear that there will be k terms corresponding 
to normal BFKL evolution where only one dipole is evolving and k — \ are 
spectators. Furthermore, there will be k{k — 1) splitting terms, according to 
the terminology we just defined, which will be proportional to the — 1-th 
density where the n—2 dipoles will be spectators. 

These splitting contributions, or to be more precise their analogue, were 
not included in the JIMWLK equation. We already start to see the conse- 
quences, since it is trivial to show that the hierarchy of equations obeyed 
by the density moments is not consistent with any sort of factorization, for 
example 



\i ^u-iVi;U2V2 j 



(46) 



This factorization is "broken" because of the splittings terms which become 



important in the region 



n ~ 



1 (or equivalently when T ~ af), that is. 



when fluctuations in the number density of particles become a significant 
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effect. We need to say here that in the original formulation of the dipole 
picture the equation for the dipole-pair density (and the higher moments) 
was not written as given in Eq. (45), but in an equivalent form in which 
the fluctuations were more difficult to recognize. Nevertheless, based on 
that picture, these low-density fluctuations and some of their consequences 
were in fact "seen" in the numerical simulations of the wavefunction of an 
evolved dipole and of the approach to unitarity^° [36, 37]. 



10. Splittings of Pomerons, Large- AT^ Equations and the 

Langevin Equation 

Now that we have derived the equation for the dipole-pair density, it is 
not hard to transform to the corresponding equation for the amplitude of 
two given external dipoles to scatter off the target. To the order of accu- 
racy it is enough to consider that dipoles scatter in the two-gluon exchange 
approximation. This elementary scattering amplitude for two dipoles (x, y) 
and (tt, v) is calculated in Appendix B and it reads (see e.g. [107, 108]) 



iv2 _ 1 
To{xy\uv)=^^-^\^ 



{x — vY{y — u 



i2n 



8 A^2 [{x-uy{y-vY 
7V2 _ 1 

= (4^j;7^Mxy\uv). (47) 

c 

For what follows, we shall set the fraction involving the color factor equal 

to one, since we have already assumed the \axgc- Nc limit in our analysis. 
To this end, the amplitude for the projectile dipole {x,y) to scatter off the 
target will be 

{T^y) = al [ Ao{xy\uv) {riuv) + {rivu) = ^^l^l{Tuv), (48) 

uv 

where in the second part, valid for -u / we have inverted the equation 
to obtain the symmetrized dipole density in terms of the amplitude for our 
later convenience. Now let us consider the scattering of a pair of dipoles 
{xi,yi) and (032,2/2) off the target. Then the extension of Eq. (48) reads 

). (49) 

Clearly, in writing the above equation, we have assumed that two dipoles 
do not interact with the same dipole. We need to say that neither a large- 
Nc nor a small coupling argument justifies this assumption. Nevertheless, 

^° In Sec. 12 we will briefly explain how unitarity comes in that picture. 
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Fig. 11. Splitting contribution to the dipole pair scattering. 



such processes will be suppressed at higher energies, as they will grow like 
a single BFKL Pomeron. Now by differentiating Eq. (49) and using the 
last, linear in (n), term in Eq. (45) for the evolution of the dipole pair 
density (the bilinear terms give rise to the normal BFKL evolution of (TT) , 
which we already know how to write), we obtain the splitting contribution 
to the evolution of the dipole-pair amplitude, which reads [38-40, 43] (see 
also [109,110]) 



dY 



split 



27r) 2tt / -^-"-"^-^oCa^iyilit^) A(a^22/2|2;f) 



xVivliT^^). (50) 



Notice that the poles of the kernel cancel with the zeros of the dipole-dipole 
scattering amplitude. The pictorial interpretation of this equation is given 
in Fig. 11. A dipole {u, v) of the target splits into two new dipoles {u, z) and 
(jz,t>), leading to the dipole kernel in Eq. (50). The two daughter dipoles 
scatter with the two projectile dipoles, as represented by the two ^o's in the 
equation. Finally the last term is proportional to the initial dipole density, 
which has been expressed in terms of the single scattering amplitude by 
making use of Eq. (48). 

This is the term that we would like to add to the r.h.s. of the second 
(large- iVc) Balitsky equation. It corresponds to the splitting of one Pomeron 
into two, and this is really the mechanism for the generation of correlations 
in higher-point functions which have important consequences to the subse- 
quent evolution of the system. For example, through Eq. (50) the single 
scattering amplitude will give rise to correlations in the dipole-pair scatter- 
ing amplitude and their significance will be realized when they feedback to 
the single amplitude through the first Balitsky equation. 

It is trivial to add the corresponding splitting term in the K-th equation 
of the large- A'^c Balitsky hierarchy. In the r.h.s. of the evolution equation for 
the amplitude for k dipoles there will be k(k — l)/2 terms proportional to 
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the amplitude for k—1 dipoles, and where in each term k—2 dipolcs will be 
simply spectators. Therefore, the new hierarchy can be easily inferred from 
Eqs. (3), (12) and (50), since in the n-th evolution equation the relevant 
processes are 1 — >^ 1, 1 — >^ 2 and 2^1 Pomeron transitions, with all other 
Pomerons being just spectators, and one just needs to take into account all 
the possible permutations. Thus, it is not necessary to present the full set of 
equations here (we shall give an equivalent compact form in Eq. (52) below), 
rather we shall indicate only its structure by suppressing the transverse 
coordinates. It reads 



dY 



= Kas (T«) -Kas (r«+^) + ^^—^ as (T^"^) , (51) 



where the K-dependent coefficients stand for the number of terms arising 
from the permutations. 

For the reasons analyzed in Sec. 9, it is clear that our hierarchy is not 
consistent with any sort of factorizing solution; the low density behavior of 
the theory has totally changed with respect to the theory without Pomeron 
splittings. Therefore, one has to find alternative ways to deal with this set 
of equations. One way to do that, is to reformulate the problem as a single 
stochastic equation. Indeed, the Langevin equation 



1 dT^y 
as dY 



2^ J -^xyz [Txz ~l" ^zy ^xy T^zTzy] 



+ £vi/ -^o(a;y|nz) Jj, V^l^lTuv v{uvz-Y), (52) 



where the noise satisfies 

{v{uiViZi;Y)u{u2V2Z2;Y')) = (5„it,2^^,iw2'^ziz2'^("s^-"s^')> (53) 

and where all other noise correlators vanish, gives an equivalent description 
[39] . We show this equivalence in Appendix I for a simple zero-dimensional 
model, while the generalization to the QCD problem at hand is straightfor- 
ward^^. However, because of the complexity of the noise correlation, this 
form may not be the best option in the search for numerical solutions. Nev- 
ertheless, one can rely on certain approximations to gain a first idea on the 
new features of the evolution. Assuming that the elementary dipole-dipole 



We should note here, that also the JIMWLK equation can bo reformulated as a 
Langevin problem [111]. However, in that case the physics is totally different since 
the noise describes color fluctuations, rather than particle number fluctuations which 
is the case here. 
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scattering amplitude is local in transverse coordinates, performing a coarse- 
graining in impact parameter space and defining the Bessel transformation 
if of the scattering amplitude one arrives at [38] 



with iy{k, Y) a Gaussian white noise, i.e. the only non-vanishing correlator is 
{u{ki,Y) u{k2,Y')) = 5{asY—asY') (5(fc^ — fc^) ^i; and where c is a constant 
of order 0{1). Notice that, up to an overall normalization factor of order 
0{l/as), ^{k) is the unintegrated gluon distribution. A numerical solution 
to this equation has been given in [112]. Furthermore, if one performs a 
saddle point approximation to the BFKL kernel in Eq. (54), the resulting 
equation is the stochastic FKPP equation which has been studied numer- 
ically in [113]. So long as the energy dependence of Qg is concerned, the 
results from these numerics are consistent with the ones from the analytical 
approach that we will present in Sec. 15. 

It is appropriate to make here a few important comments on the hierar- 
chy we have derived, and which is compactly written in the "exact" Langcvin 
Eq. (52). The first is with respect to the strength of each term. In the low 
density region where T ~ a^, one can see in the Langevin equation that 
both the noise (splitting) term and the BFKL terms are of the same order, 
while the merging term is suppressed by a factor C(a^). In the intermediate 
region where ^ T <C 1, the BFKL terms dominate, and in the region 
near the unitarity limit T ~ 1, the BFKL terms and the merging term are 
of the same order, while the noise term is subdominant by a factor C(a^)^^. 
The second comment wc should make is that, as in the case of the Balitsky 
equations, T = 1 is still a fixed point of the evolution; the noise term in 
the Langevin equation vanishes for constant T due to the presence of the 
Laplacians. We should also say that we do not expect the term we added to 
describe properly the Pomeron splittings in the high density region, since 
we have heavily relied on the two-gluon exchange approximation. In fact, 
even though we expect on general grounds the Pomeron splittings to be 
positive, the term we added can be negative in some regions [115]. Never- 
theless, this is not very worrisome, since Pomeron splittings are supposed to 
dominate in the low density region, and to this end the leading contribution 
has been taken into account. In that region the r.h.s. of Eq. (50) is clearly 
positive, since V^V^Tu^u is in fact proportional to the dipole density which 

Terms which arc of order 0{6lsT^) when inserted to the r.h.s. of Eq. (52) have also 
been calculated [114]. One can see that such term are subdominant in all regions. 




P 



+V2ca2<^(fc) u{k,Y), 



(54) 
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is positive, while away from that region the term is anyway suppressed with 
respect to the other contributions as we have just seen. 

Before closing this section, and in order to cross-check that our ap- 
proach was correct, let us compare with a well-known result obtained from 
the construction of Pomeron vertices in perturbative QCD. By neglecting 
the non-linear evolution terms, it is clear that one arrives at a solvable prob- 
lem. Thus, by first solving the BFKL equation, and subsequently solving 
the second equation in the hierarchy, which is an inhomogeneous one, and 
assuming the initial condition (TT) = 0, we find [39] 

Y 

{TxiyiTx^y^) = ^ J J ''^uvzAY~y{Xiyi\uz)AY-yiX2y2\zv) 

uvz 

xK'^liTuv). (55) 
This coincides with the Idxge-Nc triple Pomeron vertex as given in [31]. 



11. Loops of Pomerons 

Since the hierarchy presented in the previous section contains both split- 
tings and mergings of Pomerons, it is clear that Pomeron loops will be gen- 
erated in the course of evolution. For example, let us find the minimal loop, 
which will be formed after a two step evolution. This is depicted in the mid- 
dle diagram of Fig. 12. In the first step the dipolc emits a soft gluon which, 
as always, can be considered as a quark-antiquark pair and the two daughter 
dipoles "radiate" two "elementary" Pomerons. In the second step, the two 
Pomerons merge into one through the emission of another soft gluon. 

Of course one understands that a "real" Pomeron loop will be formed 
after a large number of steps, if we dress all parts of the process with normal 
BFKL evolution. Following this procedure we can construct the Pomeron 
loop vertex and let us make a short digression here, to say that one can do 
the same for all kinds of such Pomeron vertices (i.e. vertices integrated over 
a large rapidity interval) that one can imagine, e.g. for the triple Pomeron 
vertex given at the end of the previous section. Of course this is a very 
tedious procedure, and at the end one has to embed all possible allowed 
combinations of these vertices into a huge diagram, which has the two initial 
colliding objects at its two ends. Furthermore, given the non-linearities of 
the system, it is not clear whether this approach to the problem can be 
successful or not. Therefore it seems that, instead of trying to calculate 
building blocks of the theory, it is more advantageous to deal with the full 
theory as a whole, i.e. the hierarchy of equations we just presented, or its 
generalization to finite- A'^c (when the latter becomes available). Of course. 
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in order to construct the evolution equations it is necessary to know the 
"reggeized" gluon transition vertices (i.e. the vertices for a single step in 
rapidity) [9,30]. In the large- A'^c limit, and for the specific observables that 
we were interested in, we only needed to have the 2 — > 4 vertex in coordinate 
space which is simply the dipole kernel Mxyz- 

Let us now return to the calculation of the minimal Pomeron loop, since 
the result will reveal some very interesting features. In terms of the evolu- 
tion equations, and starting from the single scattering amplitude (T), the 
diagram in Fig. 12 corresponds to the generation of (TT) through the fluc- 
tuation (splitting) term given in Eq. (50), which subsequently gives feedback 
to (T) through the non-linear term of the first Balitsky equation^^. Assum- 
ing that the initial state of the target was described by a single dipole at 
{u,v) and the projectile dipole is at {x,y), it is easy to show that [43] 

— 2 r 

PL° = -2a^(^) MxyzM^^^Aoixz\uw)Aoizy\wv) (56) 

zw 

times a factor of order (AY)'^ = (step)^ which we have suppressed. We im- 
mediately see that this result is free of any possible divergencies, since the 
poles of the dipole kernels are canceled by the zeros of the ^o's- Further- 
more, and as expected, the contribution of the Pomeron loop is negative and 
leads to a decrease of the scattering amplitude(s) in the course of evolution. 

Even though we were thinking that we put the "whole" evolution in 
the wavefunction of the target dipole {u,v), it is clear that Eq. (56) serves 
for an alternative interpretation, as shown in the last diagram in Fig. 12. 
Both the original, target and projectile, dipoles split into two child dipoles 
(each of them), processes which are represented by the two dipole kernels 
times Q^. Afterwards, the ensuing dipoles scatter with each other at the 
two-gluon exchange level, and this gives rise to the two ^o's accompanied 
by a factor af. Therefore one arrives at the conclusion, that a merging of 

In terms of the language used in the forthcoming Sec. 12, this procedure corresponds 
to the successive operation -ff|_,2^2-.i on the single amplitude. 



40 



D.N. Triantafyllopoulos 



Pomerons in one of the scattering objects (here the target) can be cquiv- 
alcntly considered as a sphtting of Pomerons in the other object (here the 
projectile). Of course, one should keep in mind that whichever is the way 
we want to view the process, the soft gluon (the red horizontal gluon in 
the figures) is always emitted in the wavefunction of one of the two objects. 
This feature of exchanging splittings with mergings by jumping from the one 
wavefunction to the other, was already valid at the level of the Balitsky- 
JIMWLK evolution, where the mergings a la JIMWLK could be viewed as 
splittings a la Balitsky. But Eq. (56) contains more than that, since it allows 
for both type of processes inside the same wavefunction. Putting all these 
properties together, we shall see in the next section that we can formulate 
a precise duality condition, which is satisfied by the Hamiltonian governing 
the evolution of the hadronic system. 



12. COM Frame, Duality and a Pomeron Effective Theory 

Now we will discuss a remarkable duality property which is presumably 
satisfied by the complete Hamiltonian at high energy and at the leading 
logarithmic level. Here we shall describe the simplified case where the two 
colliding hadrons are described by non-saturated wavefunctions [43], while 
an approach to the general proof has been given in [42] , where the idea of 
the duality was explicitly written for the first time. 

So let us consider the scattering of two evolved dipoles. Assuming that 
the right (R) moving object is evolved hyY — y and the left (L) moving one 
by y, one can write the average S-matrix as [105] 



{S)= J VaRVai^ZY-y[aK]Zy[ai^]exp i j fj^{z) a'^{z) 



(57) 



The last exponential factor stands for the S'-matrix for a given state of both 
the R and L evolved dipoles, and it is just the eikonal coupling between the 
color charge in the L object and the field created by the R object. Notice 
that it can be brought into a symmetric form under the exchange R^L, 
by using the Poisson equation V\a'^i^{z) = —p^^-^{z) and integrating by 
parts. To obtain the average S'-matrix we simply integrate over all possi- 
ble configurations for both wavefunctions. This formula would be exact in 
QED, however it has a restricted range of validity in QCD. As we have seen 
in QCD, and if one wants to include all multiple scattering contributions, 
one needs to use path order exponentials to account for the color matrices 
non-commutativity in the interaction vertices, so Eq. (57) cannot be cor- 
rect in general. Nevertheless, it contains the possibility that any number 
of projectile dipoles interact with the same number of target dipoles, even 
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though each individual dipole interacts, at most, only once. This simul- 
taneous scattering of many dipoles is equivalent to a "multiple Pomeron 
exchange" and one expects the final result to respect unitarity limits. Thus, 
Eq. (57) goes much beyond the BFKL equation. We can show that all the 
above approximations and conclusions are meaningful in a frame which is 
not fully asymmetric, i.e. when y = kY, with k a constant of order 0{1), 
but of course smaller than one, for example in the center of mass (COM) 
frame y = Y/2. Requiring that the two wavcfunctions be unsaturated, one 
determines the critical rapidity Yc up to which the approach will be valid. 
In the center of mass frame one will have 



nR/L{Y/2) - al exp[a;py/2] < 1 ^ Y, 



2 1 1 



at 



(58) 



while in the case of a not fully symmetric frame one has to replace the factor 
of 2 in Yc with l/max(K, 1 — k) > 1. But the onset of unitarity will come at 
a much lower rapidity [7]. Indeed, for Y = Yc/2, with Yc given by Eq. (58), 
the contributions of the previously mentioned simultaneous scatterings is 



[a2nR(n/4)nL(V4)]''~l, 



(59) 



while in the case of a not fully symmetric frame unitarity will come at 
Y ^ max(Av, 1 — Av) ^ (l/wp) ln(l/a^) < Yc, which is of course a frame 
independent value as it should. Simply, the choice of the COM frame is the 
optimal one, if one wants to study unitarity without having to worry about 
saturation. The integer n in Eq. (59) refers to the number of Pomeron 
exchanges, for example n = 1 simply gives the BFKL contribution. By 
proper summation of all these terms, which are of the same strength and of 
order 0(1), the total S'-matrix will respect the unitarity constraints [7]. 

After this short digression to justify the use of Eq. (57), let us return to 
the proof of the duality property, which is based on two natural conditions 
[42,43]. The first is that both objects are described by the same physics, 
more precisely both wavefunctions appearing in Eq. (57) obey the same 
evolution law in rapidity. The second is that the total averaged S-matrix 
is independent of y, i.e. independent of the precise separation of the total 
rapidity interval Y, a condition which is clearly dictated by Lorentz (boost) 
invariance. Thus, by setting the derivative of Eq. (57) with respect to y 
equal to zero, and assuming that the evolution Hamiltonian has the same 
functional form for both hadrons, we easily find that 



^ ~ J ^Q'R exp 
x\zY-y[an]H 



ah, 7 



Zy[ah\-Zy[ah\H 



iJofR 



ZY-y[aK] . (60) 
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Now, by performing an integration by parts with respect to or in the second 
term, using the identity H[a,S/{i5a^)]S = H^[5/{i5pi^), pi,]S, where S is 
the exponential standing for the eikonal S'-matrix, and finaUy performing 
another integration by parts again in the second term, but with respect to 
ol this time (recall that V^ai^{z) = —pi^{z)), we find 



= Va^ Vai, exp 
xZY-y[aR] \h 



i5pi 



]zy[ai^]- (61) 



H 


' 5 ' 
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1 oa 




1 op 



The term in the square bracket must vanish and we arrive at the duality 
property 



(62) 



where in H all the functional derivatives stand to the left, while in its 
Hermitian conjugate they stand to the right. 

This duality constraint provides an alternative way to see that the 
JIMWLK Hamiltonian given in Eq. (19) can not be the end of the story. 
Indeed, while it contains an arbitrary number of color fields a" (through 
the expansion of the Wilson lines), it contains only two functional deriva- 
tives with respect to those fields. Of course Eq. (62) cannot determine the 
"complete" Hamiltonian, but it could facilitate the search for it, or at least 
help us to construct an "approximate" one which contains all the essen- 
tial physics. For example, we shall first try to find out how the evolution 
equations we presented in Sec. 10 can be put in a Hamiltonian form [43]. 

To this end, we shall recall the Hamiltonian proposed in [40] in order to 
describe the splittings of Pomerons. For our purposes it is more convenient 
to present it in the form 



H 



t 

1^2 
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(63) 



The notation 1 ^ 2 is clear, since the Hamiltonian contains two factors of 
p and four factors of 5/5p and therefore can lead to a transition from one 
to two Pomerons (and thus from two to four exchange gluons). Indeed, the 
dipole density is bilinear in the charge density, more precisely 



+ n^ 



1 



Pu Pv ' 



(64) 
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and then by acting the Hamiltonian on nuivinu2V2i which corresponds to 
the, symmetrized under the exchange of quark and antiquark legs, dipole- 
pair density (with the appropriate assumptions that the two dipoles have 
not zero size and cannot coincide), it is straightforward to show that the 
outcome is the sphtting term in Eq. (45). For this to happen, one needs to 
assume the large- A^c limit, and this means that one neglects the terms arising 
from the action of two derivatives with the same color index on two sources 
of different color. Similarly, and as a cross-check, we can act on Txiyi^Tx2y2, 
which corresponds to the simultaneous scattering of two external dipoles 
off the target. Assuming the two-gluon exchange approximation (for each 
dipole), and making use of the Poisson equation, we obtain the splitting 
term as given in Eq. (50). 

Of course this Hamiltonian is not the dual of the JIMWLK one, as 
can be seen by just comparing the relevant expressions, two of the reasons 
being that the latter does not involve any large- A?c approximation and it 
can describe any n — ^ 2 vertex, with n > 2. However, in order to build 
our approximate Hamiltonian let us construct the dual part of Eq. (63) . By 
letting p — ^ —iS/Sa and S/Sp — ^ ia we find 

uvz 

By acting on the dipole scattering amplitude Txy, again at the two-gluon ex- 
change level, we see that we obtain the non-linear term in the first Balitsky 
equation. To complete the construction within the present approximations, 
let us extract the BFKL part of the JIMWLK Hamiltonian. By expanding 
Eq. (19) to quadratic order in the color field a", integrating over the lon- 
gitudinal coordinate x~ and "taking the laxge-Nc limit by an appropriate 
contraction of color indices" [43] , we obtain 

uvz 

This part of the Hamiltonian will generate the normal BFKL evolution for 
the amplitude for a projectile dipole to scatter off the target (always at the 
two-gluon exchange level). Furthermore, it is self-dual; if we let a — > —iS/Sp 
and S/6a ^ ip, the obtained Hamiltonian will generated the normal BFKL 
evolution for the target dipole density, which is equivalent to the evolution 
of Txy as we show in Appendix C. 

By adding the three pieces given in Eqs. (63), (65) and (66), we obtain 
the total Hamiltonian 

H'' = hI + hI^^ + hI^„ (67) 
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which gives rise to the hierarcliy of equations we liavc presented in Sec. 10. 
It satisfies the self-duality condition (62) and it describes n ^ k, k^k + I 
and K+l ^ K Pomeron transitions, where, in all 1 of the Pomerons 

are simply spectators. Of course we should emphasize again that this is not 
the complete solution. Nevertheless, it contains the essential features of 
"Pomeron dynamics" and a possible numerical solution to this Hamiltonian 
problem would presumably describe the approach to unitarity in a realistic 
way. 

Before closing this section, and in order for any possible confusion to 
be avoided, we should stress that the duality transformation interchanges 
Pomeron splittings with Pomeron mergings within the same wavefunction, 
i.e. it changes the physical process. Thus, if the Hamiltonian contains both 
processes, then it will be self-dual, i.e. invariant under the transformation. 
On the other hand, a given process can be always viewed either as a Pomeron 
splitting in the wavefunction of the one hadron, or as a Pomeron merging 
in the wavefunction of the other hadron. 



Before proceeding to the quest of a complete Hamiltonian, let us try to 
construct a generalization of the splitting Hamiltonian -ff|^2 i'^^] written in 
the previous section. Recall that the JIMWLK Hamiltonian is written in 
terms of Wilson lines and can describe any n ^ 2, with n > 2, transition 
vertex at finite- A^c; as sketched in the left part of Fig. 13. Therefore it seems 
reasonable to try and find a similar way to describe all 2 ^ n vertices, as 
shown in the right part of Fig. 13, without restricting ourselves to the large- 
Nc limit. We shall not give the derivation here, but rather we shall present 
the answer as the dual of JIMWLK, which reads [41] (see also [44]) 



where the Wilson lines in the adjoint representation are defined through 



and with the sources appearing in H evaluated at light-cone time x"*" = oo. 
One way to see why the light-cone time x'^ appears in the above equations 
(while at the same time x~ is absent, or better it has been integrated over), 
is that the 2 — > n vertex inside the wavefunction of the right mover (target) 
can be also viewed as an n — >^ 2 vertex inside the wavefunction of the left 



13. The Dual of JIMWLK 




(68) 



uvz 




(69) 
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Fig. 13. Vertices of JIMWLK and its dual. 



mover (projectile). Then the appropriate longitudinal variable for this left 
mover is x^. In general when one applies the duality transformation, one 
needs to let <^ . 

Thus, given this Hamiltonian we could write the equations of motion 
for our observables. Here we shall do things in a way which is a little 
bit different than, but equivalent to, the one we followed in the JIMWLK 
case. The variables appearing in Eq. (68) are the color charges and the 
temporal Wilson lines W, which satisfy the commutation relationships 



=-i5rV.^x.., (70) 

where the first two arise from the definition of the Wilson line in Eq. (69), 
while the last is imposed dy the Jacobi identity 

[Pl.. [pI Wf] ] + [pI [Wf, pl\ ] + [Wf, [pi, pI] ] = 0. (71) 

Perhaps the non-commutativity of the color charges should not come as a 
surprise since we are interested in describing a system whose density is not 
high. In fact, close to the fluctuation line we defined earlier, the dipole 
density is of order 0{1), and therefore charge densities are of order 0{g) (of 
course {p) = due to color neutrality, but (pp) = 0{g'^) and so on). For the 
same reasons, in the JIMWLK evolution there was not such an issue, since p 
was always much bigger than and the commutator was subleading. Now 
the evolution equation for an appropriate observable O may be written in 
the form 

^-^ = {[HM)- (72) 

Based on the above considerations, it becomes obvious that if we consider 
gauge invariant observables built from the temporal Wilson lines in the fun- 
damental representation, and starting from iTiWxWy), we shall generate 
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the dual of the Bahtsky hierarchy. This is not so surprising after all, since 
the diagram on the right part of Fig. 13 when viewed upside-down corre- 
sponds to the JIMWLK evolution of a left mover, while the correlators of 
the temporal Wilson lines correspond to the amplitudes for right movers 
to scatter of these evolved left movers. But of course this trivial aspect is 
not what we were looking for, rather we are interested to see how H can 
describe the low density fluctuations. Before that, and as a first check, one 
can see that Eq. (72), with O the bilinear of the charge density as defined 
in Eq. (64), will lead to the BFKL equation for the dipole density of the 
target hadron. To measure the fluctuations we need to consider a higher 
order correlator of charges, for example the dipole pair density n^uivi ^■U2'U2 
with n defined in Eq. (64). Here, one will obtain Eq. (45) but there will be 
additional terms to the right hand side of the equation. The same will hap- 
pen if we consider the scattering of two projectile dipoles in order to probe 
these fluctuations; acting with H on Tx^y^Tx2y2 at the two-gluon exchange 
approximation, we obtain the normal BFKL evolution terms, the splitting 
term in Eq. (50) plus additional terms. Some of these terms will be sup- 
pressed in the large- A^c limit, but some will not. For example there will be 
terms which correspond to the simultaneous scattering off the same target 
dipole, and this process is not subdominant in the multi-color limit. More 
generally the two dipoles could scatter off more complicated combinations 
of the target color sources. 

One may have already seen that there are some subtleties associated with 
the dual Hamiltonian and the particular observables that we are interested 
in. Since the color charges are non-commutative, one has a difficulty in 
defining the observables at equal time x^. For example, even though there 
is no problem with the single dipole density, since [p^, p^] = 0, there is an 
ambiguity in defining the dipole-pair density since [fiuivi , nu2V2] / and 
similarly for the dipole-pair scattering amplitude, since the color field a" 
is related to the charge through the solution of the Poisson equation. 
Thus one needs a prescription to circumvent the problem, for example by 
considering symmetrized observables. 

Of course these difficulties should disappear when one assumes that the 
target be composed of dipoles. The dipoles are color neutral, and therefore 
there should be no non-commutativity issues any more. In that case, one 
may ask if the dual Hamiltonian H will reduce to -ffi^2 i^^) since 

the latter was derived on the basis of the dipole picture. But one cannot 
really answer this question, since the assumption that the relevant degrees 
of freedom are dipolar, is a condition to be imposed on the wavefunction of 
the system and is not a property of the Hamiltonian. So let us assume that 
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the wavefunction is of the form 

Mp] = E / d^NPN{{z^}■,Y)ll — tT{W^_,Wl)^{p), (73) 
Ar=i 1=1 

with the notation introduced in Sec. 9. Now in order to find the evolution 

equation of (C), it is much easier to obtain first the evolution of the wave- 
function Zy[p\-, then multiply with O and finally integrate over p. When 
we apply this procedure to the dipole-pair density we obtain the standard 
evolution equation given in Eq. (45), provided that during the calculation 
we drop terms which are suppressed in the multi-color limit [116-118]. Still, 
this description contains more than the one given by h\^2- the latter, 
only two gluons could be radiated from each color dipole in the wavefunction 
of the target, while in the former there is no such restriction and a dipole 
can radiate an arbitrary number of gluons. In particular, this means that 
two projectile dipoles are allowed to scatter off the same target dipole. 

14. Effective Action 

The next natural task should be to construct an effective theory which 
takes into account all the vertices describing m ^ n transitions, with m and 
n arbitrary (and greater or equal to 2), like the generic one appearing in 
Fig. 14. Whether this is the full set of vertices contributing at the leading 
logarithmic level in asln(l/x) or not (presumably not), is a separate issue 
and if not it could turn out that the search for the "final solution" be a 
formidable and too ambitious task. In any case, it seems reasonable to 
avoid the discussion of this problem here. 

The diagram in Fig. 14, and in the spirit of the CGC approach, corre- 
sponds to the propagation of a semi-fast gluon (the horizontal red gluon) in 
the presence of two types of background fields, and it represents the modes 
that we need to integrate as we go to higher and higher values of rapidity. 
The two background fields (in the Coulomb gauge) are the component 
of the color field, represented by the vertical gluons above the semi-fast one, 
and the radiative part of the A~ component, represented by the vertical 
gluons below the semi-fast one. The semi-fast gluon will be slow with re- 
spect to the component of the background field, so it will appear to the 
latter as moving fast to the positive x~ direction. Thus it will couple to that 
component through a Wilson line in the longitudinal direction. Similarly, it 
will be fast with respect to the A~ component, it will appear to the latter 
as moving fast in the positive a;+ direction, and thus it will couple to that 
component through a Wilson line in the temporal direction. Therefore, and 
not surprisingly, both types of Wilson lines which we encountered (sepa- 
rately) in the cases of the JIMWLK Hamiltonian and its dual will appear 
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in the final form of our effective action, and in fact these will be the only 
degrees of freedom. In our current notation these Wilson lines in the adjoint 
representation (and by dropping the tilde we have been using so far) will 
read 



V^+{x) = Fexp 
T^,-(a;) = Pexp 



OO 

ig j dx"A+(x+,x",a;)T" 

— oo 
oo 

\g j Ax+A-{x+,x-,x)T'' 



(74) 
(75) 



Again we will just present the final result for the change of the effective 
action under a rapidity step, and it can be given in the form [44] 



A5, 



eff 



Tr 



+ perm. 



(76) 



where "perm" stands for the three more terms arising from the possible 
permutations in the position of the spatial derivatives, and with the indices 
in the longitudinal and temporal Wilson lines standing for their and x~ 
dependence respectively, in accordance with their definitions given above. 
This expression for the effective action has also been obtained in [45] by 
studying the scattering of two "shock waves" as developed in [119], while 
an earlier stage of the same expression has been reached in [42]. The four 
Wilson lines are not independent, rather they satisfy the constraint 



vlw^^v^^wl = 1 



iV2 



1 



TrTy^ = l, 



(77) 



where is the Wilson loop shown in Fig. 15. This loop encloses the 
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diamond defining the "interaction" area; the color field has support in 
the strip < x~ < Xy ~ ^o^^ width expands with increasing 

rapidity, while the color field A~ has support in the strip < x"^ < Xy ~ 
x^e~^ and its width contracts with increasing rapidity. Notice that the 
value of the effective action in Eq. (76) depends only on the value of the 
fields on the boundary of this diamond-shaped area, and therefore it seems 
that the high energy problem reduces to a two-dimensional effective theory. 
Of course, at the same time, one will not be able to probe x~ and/or 
x^ dependent correlations of the system (but this was the case in both 
JIMWLK and its dual anyway), since the construction of Scs results in a 
coarse-graining along these longitudinal and temporal directions. 

This effective action is invariant under duality transformations as antic- 
ipated [44]. Indeed, if we let 

Woo - Vl, W.^ ^ Vl^, (78) 

Seff transforms to an expression which is fully equivalent to the one given in 
Eq. (76). Furthermore, the effective action will reduce to known expressions 
in limiting cases. Expanding the temporal Wilson lines to lowest non-trivial 
order, integrating over x"*", letting A~ —\6/6p and using the Poisson 
equation to relate p to A'^ = a, one obtains the JIMWLK Hamiltonian as 
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given in Eq. (19). Similarly, an expansion of the longitudinal Wilson lines 
leads to the dual of JIMWLK as given in Eq. (68). 

However, we should stress that Eq. (76) is not the complete answer to 
our problem. What we would really like to have is a Hamiltonian, that 
is, the effective action should be accompanied by well-defined commutators 
satisfied by its degrees of freedom, i.e. the Wilson lines. Then we would be 
able to find how the wavefunction of the system and the average values of 
the observables evolve with rapidity, as we did in the two limiting (JIMWLK 
and its dual) cases. For example, we have not been able to do the "simplest" 
but most natural exercise; to write the evolution equations obeyed by the 
dipole scattering amplitudes in this general case where all m — n transitions 
are possible, as this would allow us to generalize the hierarchy presented in 
Sec. 10. This problem is still "open" and therefore we shall not discuss it 
here. Let us just mention that the simple rule A~ — > —i6/6p for promoting 
an action into a Hamiltonian, and which is valid in the two limiting cases, 
will presumably not work any more. First of all, A~ is not the natural degree 
of freedom, rather the latter is given by the corresponding temporal Wilson 
line. Furthermore, a replacement like the one mentioned above would lead 
to commutators which depend on "broken" Wilson lines (lines in, say x~ , 
which run from — oo to an upper limit different than oo), a feature that 
makes the whole approach ill-defined. 

15. The Saturation Momentum Revisited 

In this last section we shall try to find the corrections in the energy 
dependence of the saturation momentum, which are induced by the splitting 
terms in the evolution equations. The result will be the one found in [34], 
even though we shall present a much simpler calculation while at the same 
time we will attempt to "improve" it. As was claimed in [34] and in Sec. 7, 
one needs to introduce an ultraviolet absorptive boundary to cut all the 
contributions from the region where the amplitude is such that T < a^, 
since they correspond to unitarity violating paths. Put it another way, a 
splitting of Pomerons in the target point of view, corresponds to a merging 
of Pomerons in the projectile point of view and therefore this boundary is 
like a saturation boundary for the projectile. The result of [34] has also 
been confirmed in [35], by identifying and studying the analogy to a similar 
problem in statistical physics. Indeed, as we saw in Sec. 10, the full large- 
Nc hierarchy can be reduced, under certain but justified for our purposes 
approximations, to a Langevin problem [38] which is similar to the sFKPP 
equation. This equation shares the same features as the, deterministic, 
FKPP equation, when the latter is supplied with a cutoff at a^. 

Therefore we would like to solve the BFKL equation with two absorp- 
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tive boundaries'^. The separation, in logarithmic units, between these two 
boundaries, the saturation Hne and the fluctuation hne, should be 

A = — — In ^ + O (const), (79) 
1 - 7r 

as shown in Fig. 6, since within A the amplitude will drop from a value 
of order C(l) to a value of order 0{a^), as determined from its leading 
power behavior in Q"^ (in Eq. (79) we have anticipated that the real part of 
the anomalous dimension will not be 7^ any more and its value 7r is to be 
determined). The constant in Eq. (79) is not really under control, since it 
can be affected, for example, by the precise value of T on the boundaries. 
Even though this constant should not affect the leading correction, since it 
is suppressed with respect to the leading term when ag 0, it will turn 
out to have an important influence on the results for small, but reasonable, 
values of the coupling constant. 

By making a change of variables, from (ln((5^/^^), y) to {z,Y), with 
z = ln{Q'^/Ql{Y)), one can write the BFKL equation as 



T = 0. (80) 



The above expression has been set equal to zero, instead of (l/a^) dT/dY, 
since we arc looking for lines of constant scattering amplitude T, that is, Y- 
independent solutions'^ for T. In Eq. (80), is the logarithmic derivative 
of the saturation momentum, which is to be determined. Now let us recall 
that the BFKL eigenfunctions are exponentials in z, where the anomalous 
dimension is in general complex, i.e. 7 = 7r + 171, with obvious notation. 
Then one can see that the unique real linear combination of eigenfunctions, 
which satisfies the boundary conditions T{z = 0) = T{z = A) = is 

7VZ 

r (X exp[— (1 — 7r) sin — . (81) 

In the above equation the sine has been obviously generated by the linear 
combination of the parts of the two degenerate eigenfunctions that involve 
the complex piece of the anomalous dimension, which is uniquely determined 
in terms of A and reads 7i = tt/A. Furthermore, in order for Eq. (81) to 

The BFKL equation, and more precisely the value of the hard Pomeron intercept, 
in the presence of momentum cutoffs has also been studied in the "remote" past 
[120, 121], but the origin of the cutoffs was much different. 

Of course, by making this assumption we will not be able to find the approach to the 
asymptotics, i.e. Y-dependent corrections for A^. Nevertheless, due to the "strong 
absorption", those corrections drop very fast, more precisely exponentially, in this 
two boundary problem. 
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satisfy Eq. (80) and since the saturation momentum is a real quantity, the 
constant Xg must be given by 

= 2^ with Im(A,) = 0. (82) 
1-7 

This last equation is the solution to our problem, since for a given value of 
7i (that is, for a given A or for a given ag), there will be a unique value of 
the real part 7r of the anomalous dimension which makes the logarithmic 
derivative of the saturation momentum real. Notice that the solution 
will not correspond to a saddle point of x(7)/(l — 7), in contrast to the 
single boundary case, cf. Eq. (31). 

When the separation of the boundaries is large, or equivalently the cou- 
pling is extremely small, i.e. A ^ 1 <C 1, we can expand the r.h.s. of 
Eq. (82) around the asymptotic anomalous dimension which is 7^ as ex- 
pected and as shown in Fig. 16. By setting the imaginary part of this ex- 
pansion equal to zero, we uniquely determine 7r and then the real part gives 
As. It is not hard to find that in this expansion the logarithmic derivative 
of the saturation momentum is given by [34] 



A. ^ X(7.) _ 7r^(l - ls)x"ils) ^ 4 00 _ 150 
as 21n2(a2) ' ln2(af)' 



(83) 



while the real part of the anomalous dimension reads 

7r = 7s + 



^ ^ (1 - 7.)X'"(7. 



3x"(7. 



- 0.372 -4^. (84) 



21n2(a2) 1^.(^2) 



It is obvious that as we let the second boundary to infinity, i.e. when ag 
0, we approach the results of the single boundary problem. The same is 
true for the scattering amplitude too; as A — ^ 00, Eq. (81) reduces to the 
scaling part (the first two factors) of Eq. (33) . Of course it is not surprising 
that the correction in A^ is negative, since there is absorption from the 
two boundaries. Moreover, we should notice that the denominator in the 
correction of A^ has again a simple interpretation; it is proportional to the 
square of the "effective" transverse phase space, in logarithmic units, for 
the evolution of the system, which here is A^ ~ \n^{a'^), as it was in the 
single boundary case (cf. Eq. (32) and the discussion after that). A rather 
important feature of these corrections, which are generated in reality by 
the negative contribution of the Pomeron loops formed in the course of the 
evolution, is that they appear to be much more significant that any next to 
leading order correction; indeed the latter would simply add a term of order 
0{as) with respect to the leading term. 
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Fig. 16. The anomalous dimension in the two boundary approximation (with the 
thin dotted hue corresponding to Eq. (84)). 

These corrections have the nice feature that they are universal, that is, 
they do not depend on the precise width of the boundary. If we let as — > kOs, 
with K = 0(1), which amounts to the modification of the boundary width 
by an additive constant, the induced correction in Eqs. (83) and (84) will 
be of order 0(1/ In^(a^)). But within the same reasoning, these equations 
are as far as one can go by using analytical methods, since any higher order 
correction will depend on the precise value of the boundary width, and thus 
one cannot rely on the BFKL equation any more. 

This brings us in a somewhat difficult situation, since the coefficient 
in the correction of is huge, mostly due to the large value of x"(7s)- 
Indeed, for reasonable values of as the correction will dominate the leading 
contribution; as given in Eq. (83) will be negative so long as Us > 0.06. 
Of course this does not mean that there is something wrong neither with 
Eq. (83) nor with the problem at hand. Eq. (83) is really supposed to be 
valid at very small a^, while for a realistic value of as we should rely on a 
(numerical) solution of the full set of equations, i.e. including non-linearities 
and fluctuations. 

One can already see that the situation will be better if we decide not 
to expand the BFKL kernel in the diffusion approximation, as we did while 
going from Eq. (82) to Eq. (83). Eq. (82) is just an algebraic equation which 
can be easily solved numerically for a given, but arbitrary, value of as- In 
Fig. 16 we show the solution for the anomalous dimension in the complex 
plane, while in Fig. 17 we show the dependence of Xs/as on as- 
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As we see in Fig. 16 and also in Eq. (84) there is an important correction 
to the real part of the anomalous dimension, but not as huge as in the case 
of Xs- It is amusing, but also natural, that there is no solution for values 
of 7 such that 7r > 7P = 1/2. This is the point where the saturation line 
becomes identical to the Pomeron intercept line. That means that the latter 
becomes a line of constant amplitude (like the saturation line); the negative 
corrections induced by the boundaries will exactly cancel the leading value 
wp = 4as In 2. Since the Pomeron intercept line corresponds to the line of 
fastest increase, we conclude that there will be no high energy growth any 
more, therefore no saturation, and thus no solution to our problem! 

In Fig. 17 we see that starts from the, well-known by now, asymptotic 
value 4.88as for vanishing as and it decreases to become zero at a point 
which is of course the same as the point in Fig. 16 where 7r = 1/2 and 
there is no solution any more. The change in the "speed" of the saturation 
momentum is still significant, but the total result is positive in a much wider 
region of possible values of a^, more precisely one has As (as < 0.3) > 0. 
Of course this is not really under control, since one is free to rescale Og 
by a factor of order 0(1). Nevertheless, this result behaves much better 
than Eq. (83), while both are solutions to the same problem and under 
the same definition for the boundary width A as a function of Og- Thus, 
one expects that in a (numerical) solution of the full problem, the outcome 
will presumably be well-defined for all reasonable values of the coupling 
constant. 
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Our presentation would be incomplete, if at this point we would not 
mention that the scaling behavior of the amplitude will not persist at very 
high values of rapidity [34] due to the presence of the fluctuations terms in 
the evolution equations [35,38]. Since the system is of stochastic nature, 
different events will lead to different profiles of the scattering amplitude 
as a function of and at a given fixed rapidity Y [35]. These profiles 
will be of the same form but shifted with respect to each other according 
to a probability density, which at a first approximation can be taken as a 
Gaussian with a diffusive radius which scales as \/ (XgY / In^ (1 / a^) . Averaging 
over all the events, one finds that [35] 

(T(Q^Y)) = f( ,MQV'Qg» ), (85) 

which clearly violates the geometrical scaling of the amplitude. Such vio- 
lations, and their consequences, for instance the breakdown of the BFKL 
approximation in the high momentum region [38], have been seen in the 
existing numerical solutions [112, 113], but the precise value of the critical 
rapidity that they set in (which could depend on the initial conditions) is 
not under real control yet. 



16. Epilogue 

In these lectures I have tried to give a simple introduction (and not 
only) to the basic theoretical developments, during the last decade or so, 
in the field of high energy (small-x) QCD. To this end, I have introduced, 
and relied mostly to, the dipole picture of high energy evolution and the 
Color Glass Condensate formulation of an energetic hadron, since, to my 
opinion, these are the simplest and perhaps most promising self-consistent 
approaches to the problem. Even though I have not followed a step by step 
procedure in the construction of the effective theory governing the dynamics 
at this high-energy limit, I have tried to describe as much as possible the 
underlying physical picture, and in the Appendices that follow I give the 
derivations of some relevant equations in order to have a more complete 
presentation. 

I feel that we have seen important progress in the field during the last 
year, since we were able to derive evolution equations containing Pomeron 
loops as one of their building blocks, which is really the mechanism that 
eventually leads to the unitarization of scattering amplitudes and the satu- 
ration of hadronic wavefunctions (even though for an accurate description 
of the latter one may need to go beyond the large- A'^,. limit). 

Of course there are still many open questions and things to be done. 
For instance, it would be natural to look for the extension of the evolutions 
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equations to finite- A?c which allow transitions among arbitrary numbers of 
Pomerons, even though the solution of such generalized equations may turn 
out to be not much different for most interesting quantities (in the same 
way that the BK and JIMWLK equations share, at some level, the same 
features). So, at the moment, we would like to have a better analytical 
insight to the current set of equations, and presumably this can be achieved 
by analyzing their stochastic nature, while at the same time more accurate 
solutions are needed. Clearly this is rather important if one aims to see the 
implications at the phenomenological level. To this purpose, and for a more 
complete analysis, one would really like to go beyond this hierarchy, which is 
good for the description of total cross sections, and extend it, if possible, to 
the case of diffr active scattering (such attempts, but restricted to BK and/or 
JIMWLK evolution have been already made [122-124]). On the theoretical 
side, I believe that it is important to investigate the relationship with the 
other approaches to the problem and the possible equivalence between them, 
so that at the end of the day one will have a unified and concrete description 
of the dynamics in high-energy QCD. 
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Appendix A 

The Dipole Kernel 

In this Appendix we will derive the expression given in Eq. (2) for the 
dipole emission kernel [5]. Let us look at the left diagram in Fig. 18, where 
a soft gluon is emitted by the quark line. With the 4-momenta of the 
quark and the gluon being = {p~ = p^/2p+ 0,p^ ^ og,P 0) and 

= {k~ = k'^ /2k~^ , = xp'^,k) respectively, the probability amplitude 
for this emission in light-cone old fashioned perturbation theory is given by 

Here, A stands for the polarization of the gluon, g i"e^ is the quark gluon 
vertex in the high energy approximation, l/k~ corresponds to the energy 
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y ^ y 

Fig. 18. The two diagrams contributing to the sphtting amphtude. 



denominator which is dominated by the gluon, and the square root factor 
arises from the usual normahzation of the color field. We have also used the 
fact that, in the light cone gauge A'^ = 0, the condition A;^e^ = implies 
= ex-k/k^ . It is straightforward to perform a Fourier transformation 
in order to obtain the probability amplitude in transverse coordinate space, 
which reads 



k 



exp [ik-{x — z)] 
(2^ 



^1 2igt- ex-{x-z) 

Similarly the probability amplitude for emission from the antiquark line of 
the dipole will be A^^(A;"^, z—y). Then the differential emission probability 
in the interval d^z and from k'^ to k^— dk~^ (we anticipate that the change 
in the probability will be positive when k^ decreases and therefore the 
longitudinal phase space "opens up") will be given by 

dP = -(27r)2^|7W^(fc+,a;-z)+7W^(A;+,z-^/)|^d2zdA;+. (A.3) 

a, A 

The factor (27r)^ arises when we transform the (square of the) wavefunction 
of the initial state in momentum space to the one in coordinate space. By 
using t'^t" = {N^ — l)/[2Nc) Nc/2 to sum over color indices in the large- 
Nc limit, and ^\^*x — to sum over polarization indices, we finally 
arrive at 

dP = ^ ^""'J^^ d'z dY. (A.4) 

2-K {x- zf{z-yY ^ ' 

In the above equation, we have introduced the rapidity increment dY . No- 
tice that this factor arises from the longitudinal phase space —dk^/k'^ = 
-dx/x = dln(l/x) = dY. 



Appendix B 

The Elementary Dipole-Dipole Scattering Amplitude 
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Here we will calculate the imaginary part of the dipolc-dipole scattering 
amplitude in the two-gluon exchange approximation. In this limit only 
tree diagrams need to be considered and one can simply use classical field 
equations. Let us assume that the dipole {u,v) creates a field a" which 
is "seen" from the dipole {x,y). In general, the amplitude for the dipole 
{x,y) to scatter off a classical field a" is determined by 



T{xy\a) 



9 

ANr 



X) 



a 



\y)\ 



(B.l) 



When this field is created by a source p"-, it can be determined by the 
solution of the Poisson equation in the Coulomb gauge, more precisely 



a^iz) 



z-wYp^lp^iw). (B.2) 



Then Eq. (B.l) becomes 
Tixy\p) 



9 



64Tr^Nr 



In 



(x 



w] 



Xy-wf 



p%w) 



(B.3) 



In the problem under consideration the source consists of the quark and the 
antiquark of the dipole {u,v), and therefore its charge is given by 

p"(w7) = 5 [d{w - w) - 6{w - v)] . (B.4) 

Substituting the expression (B.4) into Eq. (B.3) and by making use of 



Nl-1 



we finally arrive at 



iV2 

To{xy\uv) = ^ 



8 



In' 



[X 



vf{y 



u 



{x - uY{y 



(B.5) 



(B.6) 



Notice that this is invariant under the exchange of the two dipolcs, i.e. xy ^ 
uv and under the exchange of the quark and the antiquark within the same 
dipole, i.e. x ^ y and/or u ^ v. The latter will not be true, in general, 
when we relax the two-gluon exchange approximation. 

One can also integrate over the impact parameter of the scattering, 
i.e. over the relative distance h = {x+ y—u—v)/2 between the two dipoles, 
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average over the orientation of the two dipoles, and multiply by a factor of 
2 to obtain a total cross section 

(TDD(n,r2) = 27ra2r|(^l + ln^^. (B.7) 

In the above equation ri and r2 are the sizes of the two dipoles {x, y) and 
{u,v), r< = min(ri,r2) and r> = max(ri,r2), while we have also set the 
color factor equal to its Idjge-Nc value. When either of the dipole sizes 
vanishes, this total cross section vanishes due to color transparency. 



Appendix C 

The BFKL Equation: From Dipole Densities to Scattering Amplitudes 

The BFKL equation for the dipole density in the wavefunction of the 
target reads 



d{nuv) " 



dY 



[Muzviriuz) +Mzvu{nzv) - Muvz{nuv)]- (C.l) 



In the two-gluon exchange approximation, the amplitude for a dipole (x, y) 
to scatter of the target is given by a linear transformation of the dipole 
density, namely 

), (C.2) 

uv 

with To{xy\uv) the elementary dipole-dipole scattering amplitude at lowest 
order in perturbation theory, and when the color-factor is set to its large-iVc 
value, one has (cf. Eq. (B.6)) 



To{xy\uv) = In 



{x — v)^{y — 

{x — uY{y — vY _ 



(C.3) 



Differentiating Eq. (C.2) with respect to Y , using Eq. (C.l) and making 
some trivial interchanges of variables we obtain 

^1^ = ^ / M^v-[Mxy\'^z)+MM=^v)-n{xy\uv)\{ri^^).{CA) 

uvz 

Let us concentrate on the integration over z, and make a change of variables 
through the conformal transformation 

^^Z+fL, (C.5) 

C2 2; - 1 
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with the parameters ci and C2 given by 

x^^ — y"^ + u^^ — x — y + u — v /^^\ 

Cl = 7 T3i 7 v^i , C2 = , (C.6) 

[xu) — (t/f ) '■ xu — yv 

and where we have defined the complex variable z = zi + iz2, with z = 
{zi,Z2) (and similarly for x,y,u,v). This particular transformation for z 
does not change the functional form of the integrand, composed of the dipole 
kernel, the square bracket and the integration measure d^z, but simply 
performs the mapping x ^ u and y v. Therefore Eq. (C.4) becomes 



xy/ 



dY 



^ J Mxyz[Mxz\uv)+To{zy\uv)-Toixy\uv)]{nu^),{C.7) 



and now by making use of Eq. (C.2), we finally arrive at the equation 
describing the scattering amplitude of a single projectile dipole off the target 
in the BFKL approximation which reads 



dY 



Appendix D 

Eigenvalues of the BFKL Equation 

Here we shall find the eigenvalues of the BFKL equation in the simplified 

case where the scattering amplitude satisfies Txy = T{r), where r = \r\ = 
\x — y\ is the size of the dipole. The integration measure on the right hand 
side of Eq. (C.8) may be written as 

oo 

d^z = 27rrir2dridr2 y di i Jo{ir)Jo{iri)Joiir2), (D.l) 



with ri and r2 the sizes of the child dipoles {x, z) and (2, y) respectively, 
and where Jq is a Bessel function of the first kind. Then the right hand side 
of Eq. (C.8) (divided by ag) becomes 



"EFKL ' 



T{r) = f ^ ^ r^£ Jo{£r)Jo{£n)Jo 



x[r(ri)+r(r2)-r(r)], (D.2) 
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with e a positive rcgularization constant which at the end wih be set equal 
to zero. Since tlicrc is no mass scale appearing in the BFKL equation and 
the kernel combined with the integration measure has no dimension, one 
expects the eigenfunctions to be pure powers of the dipole size r. Indeed, 
with r(r) = r'^^^~'^^ one can successively integrate over ri, r2 and i in 
Eq. (D.2) to obtain 



r(7-2e)r(i+e-7) r(i 



^2(l-7)+4e_p3) 



.r(7-e)r(l+2e-7) 4^r(i+e)J' 
Taking the limit e — ^ one finds 

/Cbfkl ^ r2(i-T) = x(7) r^^'-^^ , (D.4) 

where 

x(7)=2V(l)-V'(7)-V'(l-7), (D.5) 

with ^(7) the logarithmic derivative of the F-function, i.e. 

, , , dlnr(7) _ , 

V^(7) = (D.6) 

An alternative, and perhaps easier, derivation may be given by noticing 
that the eigenvalue may be written as 

z 

where we have set y = 0, we have used the fact that the first two terms in 
the BFKL equation give the same result and we have split the virtual term 

into two parts using x^ = x-z — x-{z — x), where both parts contribute 
equally. Now we can make the change of variable z = xu to obtain 

1 K 7j^^^~^'^ + ^^^~^-2cos0 



2 





. — 2u cos 4) + v? 




where we have transformed the integration over u from 1 to 00 to an inte- 
gration from to 1, by letting u — > 1/u. Now the integration over u leads 
to Eq. (D.5). 
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Appendix E 

The JIMWLK Hamiltonian 

In this Appendix we "sketch" the derivation of the JIMWLK Hamilto- 
nian. Let us decompose the field which is associated with the wavefunc- 
tion of the energetic hadron as 

A'' = B'' + a'' + SAi", (E.l) 

where is the background field generated by the source p, represents 
the semi- fast modes with longitudinal momenta such that <C \p'^ \ ^ 
and with 6 <C 1, while SA^^ contains modes with momenta < 6A+. Since 
we want to measure expectation values of observables at scales below or 
equal to bA'^ , we need to integrate the semi-fast modes a'^. Therefore we 
define an effective action as 



A+ 

'eflf 



-i In J Va exp {iS [A,p]), (E.2) 
6A+ 

where the action S [A, p] generates the classical field equations in the pres- 
ence of a source p (cf. Eq. (15)) when 6A^^ = 0. The background field is 
determined by the solution to the Poisson equation. In the light-cone (LC) 
gauge, where only the transverse components are non-zero, i.e. B^^ = S'^^B^, 
we have 

D^F''+ = -A d+B' = p{x-,x). (E.3) 

In writing the above equation we have used that F~+ = —d~^dA~ ~ 0, which 
follows from the fact that SA~ contains modes with very small longitudinal 
momenta and therefore it varies very slowly with x~ . We can write the 
solution to Eq. (E.3) as 

B' = ^U{x-,x)dl.U\x-,x), (E.4) 

where the Wilson line is given by 

r ^~ 1 
C/t(a;",a;) = Pexp ig J dz' a" {z' , x) , (E.5) 
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and similarly for U. In the above equation a = is the background field 
in the Coulomb gauge and, with p the corresponding Coulomb gauge source, 
it satisfies 

Vla''{x-,x) = -p%x-,x). (E.6) 

Notice that is a two-dimensional pure gauge, i.e. F"^^ = 0. Now we 
expand the action 5 around Aq = 5'^^B'^ + 5'^~5A~ and to second order in 
the semi- fast modes a^, that is 

S = So + ^^a^ + \a>^G-lar (E.7) 

The coefficient of the Hnear term may be written as 

— = D^F"' = D^F^ + D+F-' + D-F+' ~ 2D+F-\ (E.8) 

where we have used the fact that F"^^ = and also the approximate equality 
D~F'^^ 2± D'^F~^. The latter arises from d'^SA" ~ which was justified 
earlier just after Eq. (E.3). In the quadratic term in Eq. (E.7) Gfj_i, is the 
propagator in the presence of the background field B^ (when calculating this 
propagator, the field SA~ can be set equal to zero to the order of accuracy). 
Now we can perform the integration over the semi-fast modes to obtain 
the (change of the) effective action. It is given by the four-dimensional 
double integral 

ASeff = - ^ / {^D+F-% G^i^ {2D+F-\. (E.9) 

xy 

Up to now we have been working in the LC gauge. However, the above 
expression is gauge invariant and it is more convenient to calculate it in the 
Coulomb gauge (we shall immediately drop the hat denoting the Coulomb 
gauge quantities so far, like in Eq. (E.6)). In this gauge the propagator 
reads 



i Ay 



e{x--y-)ul_^_{x) + eiy--x~)Uy-,-ix)\,{EAO) 



with Ay = In (1/6) representing the differential enhancement in the lon- 
gitudinal phase space and where the Wilson line is given by Eq. (E.5) 
but with the lower limit replaced by y~ (and similarly for the Wilson line 
U). Using the fact that this propagator satisfies D^Gxy = 0, it is straight- 
forward to show that the integrand in Eq. (E.9) is a total derivative with 
respect to both x~ and , so that the result of the integration comes only 
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from the "surface" terms. Furthermore, since the propagator is independent 
of the hght-cone time, we can immediately integrate over x"*" and . Using 
= —d^5A~ and defining A~[x~,x) = J dx^6A~ (x) we arrive at 

-5M;(oo)V;taM;(-oo)-5M-(-oo)y,5M-(oo)],(E.ll) 

where the Wilson lines are determined by Eq. (20). Now the evolution 
Hamiltonian can be obtained via the replacement A~ — —iS/6p. Using the 
Poisson equation V^a = —p we can express S/Sp in terms of 5 /5a, and 
recalling that the functional derivatives will act at the end-points of Wilson 
lines we can infer that 

iVl)'"' = Vx'Tjr-T, (E.12) 



(5a«(-oo) (5a^(oo)^ '^^ (5q;^(oo)' 

so that we can express the functional derivatives at a;~ = — oo in terms 
of those at x~ = oo. Then we obtain the JIMWLK Hamiltonian in its 
"standard" form 



uvz 



where the kernel K^vz is given by 

i^u-z)-{z-v) 

{u-zf{z-vf ^ ' 

When the Hamiltonian acts on the gauge-invariant observables defined in 
Eq. (21) one can let 2Kuvz -Muvz to arrive at the simplified form given 
in Eq. (19) [26]. 

Appendix F 

The Balitsky Equations 

Here we derive the first two Balitsky equations starting from the JIMWLK 
Hamiltonian. The latter is given by 



uvz 
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while the evolution equation for a generic operator O is determined by 

d{0) 



dY 



{HO). 



(F.2) 



First we would like to find the action of the JIMWLK Hamiltonian on the 
S'-matrix 



(F.3) 



which describes the scattering of a single dipole off the hadronic target. The 
action of the functional derivative on the Wilson lines is 



Sa"- 



(F.4) 
(F.5) 



Now we easily find that 



S 



Sa^ Sal '"^ 



El 



Sx-uSyvtloitH'^V^Vy), 



(F.6) 



where we have dropped terms proportional to S^v since the kernel Muvz 
vanishes when u = v, and we have anticipated that both terms in the square 
bracket in the above equation will contribute the same. Let us consider the 
first term in the square bracket of Eq. (F.l). Using = {N^ — l)/{2Nc), 
we find the first contribution to HS. 



xy 



Ar2 

first = - 



2m 



27r J 



M 



xyz Sxy. 



(F.7) 



Now consider the contribution coming from the second term in Eq. (F.l). 
In order to transform the adjoint Wilson lines to fundamental ones, we shall 
make use of 



Then it is straightforward to show that the contribution of the second term 
is the same as the one of the first, that is 



AT^ - 1 a 

J. T „ J. 

second = — 



27V2 



27r J 



(F.9) 
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The third term involves 

{VirV^'tr{tH'^Vj,Vy) =tT{V^t^VzV^t^Vy), (F.IO) 
and by using the identity 

tr(tMr5) = i tr(A)tr(i?) - ^ tr(^S), (F.ll) 

which arises from the Fierz identity 

we find that the corresponding contribution (which is also equal to the one 
coming from the fourth term) is 



third = fourth = - — / Mxyz 



2 27r 7 ■ 



(F.13) 



Combining Eqs. (F.7), (F.9) and (F.13)^^ we finally arrive at the first Bal- 
itsky equation 

d{Sxy) Oi, 



dY 



= g / Mxyz [{SxzSzy) - {Sxy)] . (F.14) 



Now wc would like to scatter two dipoles {x\,y-\) and {x2,y2) off the 
hadron. The corresponding S-matrix is given by 

Sxiyi;x2y2 — Sxiy\Sx2y2i (^-15) 

with Sxy the S'-matrix for a single dipole as determined by Eq. (F.3). In 
order to find the evolution equation for this operator we need again to act 
with the JIMWLK Hamiltonian. When both functional derivatives act on 
the same dipole, then the other dipole is just a spectator and the evolution 
equation can be trivially obtained as a result of the Leibnitz rule and the 
first Balitsky equation. We have 

d{SxiyiSx2y2) _ O^s f ,/<-, q _ Q ^ Q \ 

QY ~ 2-K I '"ifi^ '"i^^fi '^xiyi)'Jx2y2/ 

z 

+ {1^2} + 0(iV-2), (F.16) 



Notice the exact cancelation of the contributions which are subdominant at large- Afc. 
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where we have anticipated that the remaining terms will be suppressed at 
large- A^ci something which will be verified in what follows. Let us consider 
the terms which arise when each of the functional derivatives acts on a 
different dipole. For example, there will be a term coming from the ac- 
tion of 5/5a%^ on Vxi and the action of 5/5a\, on Vy^. For this particular 
contribution we have 



Using the identity 



j;^Sx,^Sy,,lv{t''VlVy,)tT{t'vlVy,). (F.17) 



tr(tM)tr(r5) = ^tr{AB) - tr(^)tr(5), 



(F.18) 



which arises from the Fierz identity (F.12) we find that the contribution of 
the first term in the square bracket in Eq. (F.l) is 



first 



4Ar3 27r j 



M 



^VlVy^VlVy,) 



HKVmMvM]. (F.19) 



Similarly, and by using Eq. (F.8), we find that the contribution of the second 
term is 



second = — 



1 as 
4N^ 27r 



M 



-— tr(y4y^Jtr(y4y^,)J, (F.20) 



while the third and the fourth term give 

1 a.. 



third = fourth 



4Ar| 27r 



Mx,y,z t^ykyzykyy^ylyy,) 



^tr(y4y^Jtr(y4yy,)]. (F.21) 



Putting Eqs. (F.19), (F.20) and (F.21) together we obtain 
1 a. 



M 



xiy2Z 



2tT{VlV,VlVy,V^Vy,) - tv{VlVy,VlVy,) 



Myyiylyy^yl) 



(F.22) 
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There arc seven more terms like the one in Eq. (F.22). Three terms are 
obtained when bjba^ acts on the first dipole (and therefore (5/(5a^ acts on 
the second dipole). They can be read from Eq. (F.22) with the replacement 
M.xxy2z ^ M.y^x2Zi—M.xxx2Zi—M.y^y2z- The last four terms are obtained 
when bjba^ acts on the second dipole (and therefore (5/(5a^ acts on the 
first dipole) and they can be simply read from the first four terms with 
the replacement 1 2. Putting everything together we finally find that 
the non-leading (in the number of colors) contribution to the evolution of 
{SxiyiSx2y2) is given by 



9{SxiyiSx2y2) 



dY 



NL ^ 2A''^ 27r ' L-" ''^lyaz + -^3/ia:2« ■M-x^x2Z •^yiy22j 

-HvlVy,vlVy,) - tr{Vy,vlVy,vl)] . (F.23) 

One sees that the terms in Eq. (F.23) are of non-dipolar structure and 
of order 0{1/N^) when compared to the terms in Eq. (F.16), as one had 
anticipated. The second Balitsky equation is obtained by identifying the 
antiquark of the first dipole with the quark of the second dipole. With the 
slight change of notation cci x, yi = X2 ^ z, y2 ^ y and z ^ w, we 
can write this second equation as 



dY 



d{SxzSzy) _ f _ c )Q \ 

2^ / ■"^'^xzw \\^xw^wz ^xzj'-'zy/ 
w 

j ■^^zyw {SxziSzwSwy Szy)) 



w 



1 etc 



2iV2 2n 



where the quadrupole operator Q is 



Qxzwy = ^ It^V^V^V^VyVlVz) - t^iV^Vy)] . (F.25) 



Appendix G 

Solving the BK equation in the High Density Region 

Let us derive the limiting form of the S-matrix for dipole-hadron scat- 
tering in the region ^^QY, ^ ^ Q^s^ with l/Q the dipole size, as de- 
termined from the solution of the BK equation, that is from the factorized 
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form of Eq. (26). Since we are looking for the sohition in a regime where 
the S'-matrix approaches its black-disk limit, i.e. 5 ^ 0, we can neglect 
the quadratic in S term in the BK equation. In order to do this properly, 
we need to restrict the region of integration in the transverse coordinates 
according to 

« (a; - zf, {z - yf « l/Q\ (G.l) 

The lower limit arises from the fact that is the boundary determining the 
transition from a region where the scattering is weak to a region where it 
becomes strong, while the upper limit determines the dominant logarithmic 
contribution to the r.h.s. of the BK equation. Thus, we easily find that 

^ = ¥ = -«.mq2/o^). (g.2) 

Using the leading behavior for the energy dependence of Qg given in Eq. (32), 
we can write the derivative with respect to Y as 

d_ ^ dM^l/Q') ^ _ x(7s) ,^ 

dY dY 51n(Q2/Q2) 1-^, ain(g2/Q2)- ^^-^^ 

Then it is trivial to solve Eq. (G.2) to obtain [37, 83, 90] 



{S) exp 



Mis) Q\ 



(G.4) 



It is very likely that this form remains the same when one takes into account 
the effects of fluctuations, but with a coefficient in the exponent which is 
reduced by a factor of 2 [125]. 

Appendix H 

The Diffusion to the Ultraviolet 



In this Appendix we shall show the significance of the high momenta 
contributions to the evolution. Since at these high momenta the system is 
dilute, we shall deal only with the BFKL equation and for simplicity we will 
even suppress the effects of the absorptive IR boundary. So let us consider 
the amplitude for a dipole of size 1/Q to scatter off a dipole of size l/fj,. 
The general (rotationally symmetric and integrated over impact parameter) 
solution is given by Eq. (6) with the replacement r — > l/Q. Now we would 
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like to find the solution for agY 1 and close to the momentum Qq defined 

by 



IJL^ exp 



1 



■Y 



Is 



(H.l) 



Notice that this is a line parallel to the saturation one if we neglect the 
prefactors. In the diffusion approximation, we can write the solution to the 
BFKL equation as 



nQ^Y) 



27ra^ _ , 
11^ 



1 



exp 



DsY 



(H.2) 



Now let us try to obtain the same solution, by performing two global evo- 
lution steps, from to Yi and then from Y\ to Y . Then one can write the 
solution, which we denote as T, in the form 



fiQ,Y)=J^TiQ,,Y,] 
^ 



J ^exp[a,x(7)(>^-n)-(l-7)ln(gVQ?)]' (H-3) 



which has a simple interpretation; the amplitude at Yi and for a given 
momentum Qi is evolved to Y, and then an integration is performed over the 
initial condition (at Yi). Performing the integration over 7 in the diffusion 
approximation and using the solution given above in Eq. (H.2) with Q ^ Qi 
and F — Yi, we obtain 



fiQ,Y) 



1 



1 



V^fDTn ^TrDs{Y-Yi) 



00 

xy^exp 




In^(QVQf) In^Qj/Ql) 
Ds{Y-Yi) 



DsYi 



(H.4) 



It is straightforward to show that the integration over Qi will lead to a 
result identical to the one in Eq. (H.2). But say we want to find the region 
in Qi, that contributes the most in order to get Eq. (H.2). Imposing an 
ultraviolet cutoff at Quv and setting, for example, Yi = Y/2, we find 



T(Q,Y) 
T{Q,Y) 



1 erfc 

2 



2ln{Ql^/Ql)-ln{QyQl) 



(H.5) 
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where erfc(a;) is the comphmentary error function, for which we recall that 



erfc(x) = < 



2 - 
1 

exp( 



exp(— x^) 
■s/ttx 



for X <C — 1 
for X = 

for X S> 1. 



(H.6) 



If we want to calculate the amplitude, say close to Qq, 
the argument of the error function can be neglected. 



results T{Q,Y) and T{Q,Y) will agree so long as ln((5uv/'3o 



the second term in 
Therefore, the two 



Thus, even if at the final rapidity Y we are interested in the amplitude 
very close to the "central line" Qq, at the intermediate steps of evolution 
all values of Q that extend up to the diffusion radius contribute to the final 
result. 



Appendix I 

The Langevin Equation and the Hierarchy 

Let us show in a simplified zero-dimensional case, how a particular 
Langevin equation is equivalent to an infinite hierarchy describing splitting 
and merging processes. Consider the equation 

an- f3n^ + = A + Bu, (I.l) 



dY 

where iy{Y) is a Gaussian white noise; {ly^Y)) = and {i^iY) I'iY')) = 
6{Y — Y'), and with all other higher noise correlators vanishing. This 
Langevin equation should be understood with the Ito prescription for the 
discretization of time, namely, if one writes Y = je, where j is a non- negative 
integer and e the time step, Eq. (I.l) should read 

^^^±^ = A, + B,u,^,, (1.2) 

with (I'j) = and (fjj^j) = (l/e)(5jj. Then it is a simple exercise to show 
that the evolution of the expectation value of an arbitrary function F{n) is 
determined by the equation 

^^ = {AF'in)) + l{B'F"in)). (1.3) 

With F{n) = n'^ and with A and B as defined in Eq. (I.l), the above 
equation leads to 



dY 



= a K (n^) - f3 K {n^+') + J \ (n^~')- (L4) 
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The generalization to the, two-dimensional and with non-local vertices, 
large-iVc QCD case is straightforward. 
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